AD-A183  £32 


UNCLASSIFIED 


A  CRITICAL  REVIEW  OF  THE  DEVELOPHENT  OF  SEVERAL 
VISCOPLASTIC  CONSTITUTIVE  THEORIES(U)  NAVAL  RESEARCH 
LAB  WASHINGTON  DC  J  EFTIS  ET  AL.  13  SEP  87  NRL-HR  3942 


»)*. ««'.<»*  -fa*  I.»  «,■  «.»  I.M.M.1  |  >  I*.  iViV 

♦I 


OTIC  TILE  CUPJ 

Naval  Research  Laboratory 

Washington,  DC  2037S-5000 


’l>  ’|i  ■ 


»i.t '»»  -.a  i.|  t 


NRL  Memorandum  Report  5942 


A  Critical  Review  of  the  Development  of 
Several  Viscoplastic  Constitutive  Theories 

John  Eftis*  and  Mitchell  I.  Jolles 


DTIC 

ELECTE 

OCT  2  3 1987 


Mechanics  of  Materials  Branch 
Material  Science  and  Technology  Division 

*Geo-Centers,  Inc. 

Suitland,  MD  20746 


I 


September  15,  1987 


Approved  for  public  release;  distribution  unlimited. 

87  10  14  291 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  1  ^  ^ 

REPORT  DOCUMENTATION  PAGE 

1*.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

lb  RESTRICTIVE  MARKINGS 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

2b  DECLASSIFICATION  /DOWNGRADING  SCHEDULE 

Approved  for  public  release;  distribution 
unlimited. 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NHL  Memorandum  Report  5942 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6*.  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Research  Laboratory 

6b.  OFFICE  SYMBOL 
(If  applicable) 
Code  6382 

7a  NAME  OF  MONITORING  ORGANIZATION 

6c  AOORESS  (Gty,  Stan,  and  ZIP  Code) 

Washington,  DC  20375-5000 

7b.  ADDRESS  (Gty,  State,  and  ZIP  Code) 

8*.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 

8b  OFFICE  SYMBOL 
(If  applicable) 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

8c  ADDRESS  (Gty,  Sun,  and  ZIP  Coda) 
Arlington,  VA  22217 


1 1 .  TITLE  (Includa  Security  Classification) 

A  Critical  Review  of  the  Development  of  Several  Viscoplastic  Constitutive  Theories 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO. 

61153N22 


PROJECT 

TASK 

NO  RR022- 

NO 

01-48 

12.  PERSONAL  AUTHORS) 

Eftis,  John  and  Jolles,  Mitchell  I. 


13a.  TYPE  OF  REPORT  ll  3b.  TIME  COVERED 

I  FROM _  TO 


16.  SUPPLEMENTARY  NOTATION 


14.  DATE  OF  REPORT  (Year.  Month,  Oay)  IlS  PAGE  COUNT 

1987  September  15  |  81 


COSATI  CODES  I  18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

GROUP  I  SUB-GROUP  |  ^Viscoplasticity  ,  s \  . 

'^'Constitutive  relations  — - -  ^ 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

‘  The  need  to  model  the  rate  dependence  that  is  inherent  in  non-recoverable  deformation  of  structural  materials 
and  the  continued  progress  in  computational  capability  have  served  as  stimuli  for  the  rapid  development  of  viscoplas¬ 
tic  constitutive  theory.  The  structure  and  development  of  these  theories  can  be  separated  into  two  categories:  those 
with  or  without  a  yield  condition.  All  theories  reviewed  here  have  common  underlying  features  and  limitations: 
phenomenological,  homogeneous  material  continuum*  unrecoverable  deformation  is  a  constant  volume  process, 
deformation  may  be  decomposed  into  linearly  additive  elastic  and  inelastic  parts.  The  theories  of  Bingham. 
Hohenemser-Prager,  Perzyna,  Chaboche,  Bodner-Partom,  Walker,  Kremple-Liu,  Miller  and  Valanis  are  reviewed. 
Comparisons  are  made  among  predictions  of  a  modified  Chaboche  and  Bodner  theories,  and  experiments  of  Inconel  at 
1200?F  for  various  loading  conditions.  /  J' 


20  DISTRIBUTION !  AVAILABILITY  OF  ABSTRACT 
CjjUNCLASSIFISO/UNLIMITED  □  SAME  AS  RPT 


22*  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Mitchell  I.  Jolles 


00  FORM  1473, 84  mar 


ACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

AS  RPT  □  OTIC  USERS  UNCLASSIFIED 


22b  TELEPHONE  (Indude  Are*  Code)  22c.  OFFICE  SYMBOL 
(202)  767-6451  Code  6382 


83  APR  edition  may  b*  us*d  until  *>h«usted  SECURITY  CLASSIFICATION  OF  THIS 

All  other  editions  are  obsolete  — 


TABLE  OF  CONTENTS 


INTRODUCTION 


(i)  Background  Considerations . 

(ii)  Index  of  Frequently  Used  Symbols  and  Definitions  . 

PART  I — RIGID/VISCOPLASTIC  AND 
ELASTIC/ELASTIC-VISCOPLASTIC  THEORIES 

Definition  of  Yield  . 

Bingham  Theory  for  Rigid- Viscoplastic  Simple  Shear  . 

Generalization  by  Hohenemser  and  Prager  . 

Perzyna  Theory  . 

4a^  Inelastic  Incompressibility  and  Isotropic  Hardening^ . 

4b.  ^Uniaxial  Constitutive  Equations  j . . 

Chaboche  Theory  . ...." . 


Sa.  xAJniaxial  Constitutive  Equations; . 

5b. '-^Modification  for  the  StrainHate  Dependence 
of  Initial  Yield  . ,ui. . 


PART  H-VISCOPLASTIC  THEORIES  WITHOUT  A  YIELD  CONDITION 

Bodnor-Partom  Theory  . 

6a.  Wodifiation  for  Directional  Hardening:  Stouffer-Bodner  Theory  . 

6b.  Uniaxial  Constitutive  Equations  ...y/.. . 

Walker  Theory . y.". . . 


7a.  Integral  Theory  . . . «. . 

7b.  Rate  Theory  . 

8.  Kremple-Liu  Theory  . 

9.  Miller  Theory  . 

10.  Valanis  Theory  . 

10a.  Intrinsic  Time  Scale . 

10b.  General  Theory  . 

10c.  Explicit  Constitutive  Equations  for  Small  Deformations  . 

(  PART  m— COMPARISONS  BETWEEN  THE  MODIFIED  CHABOCHE  AND  THE 

\  BODNER  THEORIES  FOR  A  HIGH  TEMPERATURE  APPLICATION 


12.  Viscoplastic  Material  Parameters  for  Inconel  at  1200°F 

13.  Experimental  and  Qualitative  Comparisons  . .1 

13a.  Experimental  Comparisons  . . 

13b.  Strain  Rate  Effects  . 

13c.  Strain  Rate  History  Effect  . 


13d.  Creep  Behavior  . 

13e.  Relaxation  Behavior  . 

13f.  Load-Unload  Reload  Behavior 
13g.  Cyclic  Load  Behavior  . 


CONCLUSION . 

REFERENCES  . 

APPENDIX— Stress,  Deviator  Stress  and  Strain  Matrices  for  Uniaxial  Loading 


Accesion  For  ^  ] 

NTIS 

CRA&! 

DTIC 

TAB 

n 

Unannounced 

□ 

Justific 

j  t'on 

By 

Ds:t  ib 

:t  / 

i 

A 

^ _ J 

Aviii  /  of 
ipvrCui 


_ ft 


A  CRITICAL  REVIEW  OF  THE  DEVELOPMENT  OF 
SEVERAL  VISCOPLASTIC  CONSTITUTIVE  THEORIES 


INTRODUCTION 


(0  Background  Considerations 

The  need  to  model  the  rate  dependence  that  is  inherent  in  non-recoverable  deformation  of  struc¬ 
tural  materials,  particularly  at  elevated  temperatures,  and  the  continued  progress  in  computational 
speed  and  capability,  have  served  as  stimuli  for  rapid  development  of  viscoplastic  constitutive  theory. 
Since  the  time  when  Bingham  and  Green,  in  1919,  first  introduced  a  theory  for  viscoplastic  simple 
shear,  more  than  a  dozen  general  theories  for  viscoplastic  material  response  have  appeared. 

The  structure  and  development  of  these  theories  can  be  separated  into  two  categories.  In  one  it 
is  assumed  that  for  strong  coherent  solids  (e.g.  structural  metals)  at  low  levels  of  stress  or  small  mag¬ 
nitudes  of  strain,  the  inelastic  part  of  any  deformation  is  negligibly  small.  The  deformation  is  essen¬ 
tially  recoverable  and  rate  independent,  i.e.,  elastic.  Inelastic  deformation,  here  taken  to  mean  defor¬ 
mation  that  is  non-recoverable  (plastic)  and  rate  dependent,  begins  only  after  the  stress  state  reaches  a 
level  sufficient  to  cause  the  onset  of  yield.  An  idealized  yield  condition  is  postulated,  with  the  yield 
function  appearing  as  an  integral  part  of  the  statement  of  the  constitutive  relation. 

In  the  other  category,  inelastic  deformation  is  assumed  to  accompany  the  elastic  deformation  at 
any  level  of  stress,  however  small.  There  is  no  yield  state  of  stress  that  separates  the  occurrance  of 
the  inelastic  part  of  the  total  deformation  from  the  part  that  is  elastic.  A  yield  condition  and  a  yield 
function,  therefore,  are  not  required  in  the  statement  of  the  constitutive  relation,  affording  thereby  a 
seeming  measure  of  simplicity,  comparatively  speaking.  However,  having  dispensed  with  the  need 
for  a  yield  function,  and  with  it  the  useful  geometrical  construct  of  a  yield  surface  in  stress  space,  this 
category  of  theories  must  resort  to  different  means  for  describing  the  hardening  that  accompanies  ine¬ 
lastic  deformation.  Internal  variables  are  accordingly  introduced  for  this  purpose,  which  are  loosely 
identified  with  microstructural  mechanisms  associated  with  non-recoverable  deformation  in  poly¬ 
crystalline  materials. 

All  of  the  theories  here  considered  have  other  common  underlying  features  and  limitations  that 
may  be  enumerated  as  follows: 

(a)  All  are  phenomenological  theories,  that  is,  the  observed  mechanical  behavior  of  a  material  is 
assumed  to  be  describable  in  terms  of  macroscopically  defined  measurable  variables.  Although 
several  theories  from  the  non-yield  condition  category  do  introduce  terminology  borrowed  from 
microstructural  models  of  polycrystalline  plastic  deformation,  the  intention  nevertheless  is  to 
describe  macroscopically  observable  features  of  inelastic  material  behavior. 

(b)  All  assume  a  homogeneous  material  continuum.  Materials  with  existing  microstructural  inhomo¬ 
geneities  are  viewed,  macroscopically  speaking,  as  being  statistically  homogeneous. 

(c)  Non-recoverable  (plastic)  deformation  is  considered  not  to  be  significantly  affected  by  super¬ 
posed  pressure  of  ordinary  engineering  magnitudes  (less  than  10s  p.s.i.).  Inelastic  deformation 
is  thus  taken  to  be  volume  perserving. 

(d)  The  theories  discussed  are,  at  this  stage  of  development,  presented  in  a  form  appropriate  for 
small  order  deformation.  The  recoverable  and  non-recoverable  components  of  any  deformation 
is  assumed  to  be  decomposable  into  linearly  additive  parts. 
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(if)  Index  of  Frequently  Used  Symbols  and  Definitions 

Since  the  theories  to  be  discussed  in  their  present  form  are  confined  to  small  small  order  defor¬ 
mation,  the  current  configuration  and  the  reference  configuration  of  a  deforming  body  may  be  treated 
as  one,  i.e.,  there  is  no  need  for  distinction  between  the  two.  Also,  all  time  derivatives  reduce  to  the 
ordinary  time  derivative,  denoted  by  a  superposed  dot. 

Symbols  for  the  variables  appearing  most  frequently  in  the  discussion  to  follow  are  listed  below. 
Any  others  that  subsequently  appear  will  be  defined  as  they  appear. 

Symbol  Name 

T,  T jk  Cauchy  stress  tensor. 

o  o 

T,T jk  Stress  deviator  tensor. 

£,  Zjk  Infinitesimal  strain  tensor. 

O  O 

E,  Ljlc  Strain  deviator  tensor. 

E '  Elastic  or  recoverable  part  of  the  total  strain. 

£"  Inelastic  or  non-recoverable  part  of  the  total  strain. 

D,  D  Rate  of  deformation  tensor. 

E  Infinitesimal  strain  rate  tensor  =*  dE/dt  a  D. 

Y,  Y £  Kinematic  hardening  tensor  variable  (stress). 

o  o 

Y,  Y jk  Kinematic  hardening  deviator  tensor, 

u,  a*  Displacement  vector. 

▼  a  u,  vt  a  uk  Velocity  vector. 

/T,  II T,  IIlj  Principal  invariants  of  T. 

R  Isotropic  hardening  scalar  variable  (stress). 

P,  p  Accumulated  inelastic  strain. 

W"  Inelastic  work. 

E,  n,  ¥  Young’s  elastic  modulus,  elastic  shear  modulus,  and 

Poisson’s  ratio,  respectively. 


0 


Temperature. 


Referring  to  a  rectangular  coordinate  system,  the  trace  of  a  second-order  tensor, 
product  of  tensors  with  tensors  and  vectors  are  defined  as  follows: 

TrT  *  T**  *  Tm  +T22  +  T33 

T  •  v  -  a:  T Jkvk  -  a, 

T-D  -  C:  T jkDu  -  Cjt 

tr  (T-  D):  TjkDkj  =  Cu  . 

The  identity  tensor  1  has  the  property. 

T- 1  *  I  T  =  T,  v  l  =  l  v  =  v. 

The  transpose  and  inverse  of  the  second-order  tensor  T  are  denoted  by  Tr  and  T“l,  while 

O 

tor  tensor  T  associated  with  T  is  defined  by  the  relation 

T  =  T  -  j(rrT)l. 

The  magnitude  of  a  vector  and  a  second-order  tensor  are  specified  by 

|  v  l  -  'Jy^y  «  Vvt  yk  ,  |T|  =  ^tr  (T  Tr)  =  \JlJk  Tjk . 

The  kinematic  tensors  £  and  D  are  given  by 

E  -  j  [gradu  +  (gradu)r],  JV*  -  y  [uJJc  +  uk j  ] 

D  «j[  grad  v  +  (  grad  v)T],  Dp  -  y  [v>Jk  +  vkJ]. 

O 

The  invariants  of  the  tensor  T  and  of  its  associated  deviator  tensor  T  are  determined  by 
/T  -  tr  T 

7/t  -  y  {(trT)2  -  trT2!,  trT2  -  tr(T-  T) 

IIIr  -  detT  -4-  (2trT3  -  3tr  T2  tr  T  +  (tr  T)3 } 

6 


/.  -  0 

T 

//.  -  -  4-rr T2,  or  II'.  «  -  VI.  -  t/T2  -  tr(TT) 
T  2  T  T 

III.  -  4-tfT3.  or  III'.  -3///.  -trT3  -  rr(TTT). 

T  3  t  T 


The  following  derivative  formula  appears  frequently  in  plastic  theory 

dir.  air. 


the  inner 


the  devia- 


(ii.10) 


(ii.ll) 


PART  I 

RIGID/VISCOPLASTIC  AND  ELASTIC/ELASTIC-VISCOPLASTIC  THEORIES 


1.  Definition  of  Yield 

On  the  basis  of  what  is  always  observed  from  the  stress-strain  data  obtained  from  uniaxial 
tension-compression  tests  of  the  structural  metals,  at  some  level  of  stress  (or  possibly  for  some  small 
range  of  stress  values)  the  deformational  behavior  of  the  material  changes  from  essentially  recoverable 
behavior  to  combined  recoverable  and  non-recoverable  behavior.  This  fact,  observed  for  uniaxial 
stress,  can  be  generalized  to  general  conditions  of  stress.  Thus,  the  state  of  stress  T  at  any  point  of  a 
deformed  body  at  which  recoverable  (elastic)  material  response  passes  into  partially  non-recoverable 
(elastic-plastic)  response  defines  a  yield  condition  at  that  point,  which  can  be  expressed  mathemati¬ 
cally  by  a  relation  of  the  form 


P 


/CD  =*• 


(l.D 


where  k  is  a  material  parameter.  Equation  (1.1)  represents  an  initial  yield  condition.  The  locus  of 
all  such  possible  states  of  yield  stress  at  any  point  of  the  deformed  body  forms,  or  defines,  an  initial 
yield  surface  in  stress  space 


F(J,k)  =  0, 


(1.2) 


that  encloses  the  origin  (zero  stress). 


For  all  strong  coherent  solids  (e.g.,  crystalline  and  polycrystalline  solids)  it  is  an  observed  fact 
that  during  loading,  increments  of  strain  are  accompanied  by  non-negative  increments  of  stress 
(presuming  no  material  micros tructural  disintegration  over  this  stage  of  the  loading).  This  charac¬ 
teristic  is  referred  to  as  stable  material  behavior.  For  such  behavior  it  is  possible  to  prove  analyti¬ 
cally  that  the  yield  surface  must  always  be  a  convex  surface  [1,2].  This  characteristic  of  a  yield  sur¬ 
face  has  been  demonstrated  experimentally  many  times  over  the  past  thirty  years. 

When  the  material  response  is  isotropic  up  to  the  point  of  initial  yield  the  yield  function  /  must 
be  an  isotropic  scalar  valued  function  of  the  invariants  of  the  stress  tensor. 


/(/T,//T,///T)  =*. 


(1.3) 


For  low  to  moderate  levels  of  pressure  (e.g.,  less  than  105  psi)  structural  metals  exhibit  very  little 
pressure  sensitivity  with  respect  to  yield.  Consequently,  as  a  first  approximation,  it  is  assumed  in 
plasticity  theory  that  the  mean  pressure 


—  / 

3  T 


(1.4) 


has  no  influence  upon  inelastic  deformation.  Equivalently,  inelastic  deformation  is  taken  to  be 
incompressible  or  volume  preserving,  [cf.  comments  at  the  end  of  this  section].  Since 


!T  +-7- (rr  T)1 


T-pl 


(1.5) 


neglect  of  the  mean  pressure  allows  the  invariants  of  the  stress  tensor  to  be  replaced  by  the  invariants 

of  its  deviator  tensor.  Furthermore  since  L  =  0,  plastic  incompressibility  reduces  the  initial  yield 

T 

condition  to  the  form 

flLJIL)  =  k.  (1.6) 

T  T 

Experimental  observations  of  initial  yield  obtained  from  combined  tension-torsion  tests  have  shown 
that  the  following  von  Mises  form  for  the  yield  function  appears  to  offer  the  best  correlation  of  the 
test  data. 


f(II'o,k)  =  JWo  -  ^2 k 
T  >  T 


where  k  =  (1/^3)  Y,  and  Yt  is  the  yield  stress  in  uniaxial  tension,  or  where  k  —  Ys  which  is  the  yield 
stress  in  simple  shear.  Note  that  the  third  invariant  of  the  deviator  stress  is  assumed  to  have  a  negli¬ 
gible  effect  on  the  initial  yield.  The  von  Mises  yield  condition  states  that  initial  yield  begins  when  the 
elastic  distortional  strain  energy  density  reaches  a  particular  value. 

The  above  expression  (1.6),  or  the  particular  form  (1.7),  for  yield  are  appropriate  for  the  com¬ 
mencement  of  yield  only.  Thereafter  yield  is  accompanied  by  hardening  which,  in  general,  is  aniso¬ 
tropic.  The  simplest  reasonably  effective  description  of  hardening  which  includes  the  Bauchinger 
effect  with  load  reversal,  is  a  combined  isotropic  and  kinematic  hardening.  For  the  isotropic  part  the 
constant  k  is  replaced  by  a  scalar  function  R  of  the  inelastic  work,  W",  or  of  the  accumulated  inelas¬ 
tic  strain,  P ,  [3,4].  For  the  anisotropic  component  a  kinematic  form  of  hardening  is  most  often 

o  o  o 

employed  which  introduces  the  tensor  Y  as  an  argument  for  the  yield  function  F  in  the  form  T  -  Y 

O 

[5,6].  The  tensor  Y  represents  the  center  of  the  migrating  yield  surface,  while  the  scalar  variable  R 
is  a  measure  of  its  size.  Both  hardening  variables  have  the  dimension  of  stress. 

Thus  for  in  the  discussion  of  yield  the  known  effects  of  the  strain  rate  and  of  the  temperature  at 
which  yield  takes  place  have  not  been  considered.  Both  effects  are  non-directional  and  can  therefore 
be  introduced  through  the  isotropic  hardening  variable  [7-11].  Keeping  in  mind  the  assumed  plastic 
incompressibility  assumption,  the  form  of  the  yield  function  appearing  in  Eq.  (1.2),  if  it  is  to  account 
for  combined  isotropic  and  kinematic  hardening,  strain  rate  and  temperature,  must  have  the  form 


where 


F(T,Y,R)  -  f(T  -Y)  -  R  , 


R  =  R{W" ,  E,0),  W"  =  fo  fr(TJE") 


R(P,±,8)  ,  P  d  |E"  |  . 


(1.10) 


The  assumption  of  the  incompressibility  of  plastic  deformation  is,  as  indicated  at  the  outset,  an 
approximation.  The  degree  of  validity  of  the  approximation  varies  with  the  structure  of  the  material 
and  the  level  of  the  mean  stress  imposed.  Experimental  data  can  be  found  in  references  [12-15], 
while  discussions  of  the  theoritical  implications  are  given  in  [16-18]. 


2.  Bingham  Theory  for  Rigid-Viscoplastic  Simple  Shear 


Bingham  introduced  the  first  theory  of  what  can  be  interpreted  as  viscoplastic  behavior  for  a 
solid  in  1919,  [19,20],  when  he  proposed  a  model  for  the  simple  shear  of  a  material  that  differs  from 
a  viscous  fluid  in  that  it  can  sustain  a  shear  stress  when  not  deforming,  i.e.,  is  at  rest. 

The  theory  assumes  that  the  material  is  rigid  and  begins  to  deform  in  simple  shear  non- 
recoverably,  when  an  imposed  shear  stress  exceeds  a  critical  value,  Ys  such  that  the  rate  of  shear 
deformation  is  proportional  to  the  excess  of  the  shear  stress  over  the  yield  stress.  Thus 

-  C[  I rI2  I  -  YSJ  for  \Tl2\>Ys, 

Dx  2=0  for  \Ta\*Yt.  (2.1) 


The  material  remains  rigid  until  the  imposed  shear  stress  exceeds  the  material’s  inherent  shear  yield 
strength  Ys .  Alternatively,  introducing  the  yield  function 


F  =  1  - 


(2.2) 


the  constitutive  relation  (2.1)  for  simple  shear  can  be  re-expressed  as 


12 


<F>  = 


IF  for /r>0,  °r\Tx2\>Ys 
I  0  for  F  s0,  or  |  7*l2 1  s  T, , 


(2.3) 

(2.4) 


with  tf  interperted  as  a  viscosity  coefficient.  The  material  parameters  Ys  and  are  temperature 
dependent.  Rewriting  the  above  relations  in  the  form 

|  r12 1  -  Ys  +  2tj£)j2,  F>0  or  |  Tn 1  >  Ys ,  (2.5) 


we  see  that  in  the  absence  of  a  shear  yield  stress  (the  material  is  no  longer  a  solid)  Eq.  (2.5)  becomes 
the  constitutive  equation  for  a  Newtonian  viscous  fluid,  whereas  if  the  viscosity  coefficient  is  zero,  it 
resembles  the  Saint  Venant-Levy  constitutive  equation  for  a  rigid  perfectly-plastic  solid  subject  to  sim¬ 
ple  shear. 

3.  Generalization  by  Hohenemser  and  Prager 


The  Bingham  theory  for  rigid  viscoplastic  simple  shear  was  generalized  to  arbitrary  states  of 
stress  by  Hohenemser  and  Prager  in  1932  by  modifying  the  yield  function  and  constitutive  expressions 
(2.2)  and  (2.3)  to  the  forms 


F 


(3.1) 


where  it  has  the  same  interpretation  as  given  in  Eq.  (1.7),  and 


1  ° 
D"  --i-  <F>  T 

2i? 


<F>  = 


F  for  F  >  0 
0  for  Fs  0. 


(3.2) 

(3.3) 


where  k  and  ij  are  temperature  dependent  material  parameters.  The  “flow”  equation  (3.2)  implies 
plastic  incompressibility  since 

tr(2rf  D”)  =tr(F  T)  =  F  tr  T  =  0  -  tr  D”  =  divv"  =  0 


which  is  the  condition  for  constant  volume  deformation.  Eq.  (3.2)  together  with  (3.1)  can  be  used  to 
obtain  the  generalized  equivalent  of  Eq.  (2.5). 


O 

T 


D"  +  2ij  D", 

v  T 


F  >0 


T  =  0  ,  Fs  0, 


(3.4) 


Coaunents: 

1.  The  form  of  the  yield  function  (3.1)  implies  no  hardening.  The  theory,  therefore,  is  for  rigid 
perfecdy-viscoplastic  flow  which  commences  when  the  value  of  F  becomes  positive. 

2.  When  if  »  0  Eq.  (3.4),  duplicates  the  Saint  Venant-Levy  equation  for  rigid  perfectly-plastic 
flow,  with  the  exception  that  in  the  perfectly  plastic  theory  flow  commences  when  F  =  0,  i.e., 
for  yj II’ o  =^2k,  which  is  the  von  Mises  yield  condition  (1.7).  For  the  viscoplastic  theory  the 

onset  of  flow  requires  that  F  >0  or  yj II’ o  >  '^2k.  In  the  absence  of  a  yield  stress,  i.e., 

v  t 

k  =  0,  Eq.  (3.4)  reduces  to  the  constitutive  equation  for  a  Newtonian  viscous  fluid. 


3.  Since  the  kinematic  variable  appearing  in  the  theory  is  the  rate  of  deformation  D,  equations  (3.2) 
and  (3.4)  are  appropriate  for  deformations  of  any  magnitude. 


4.  Perzyna  Theory 


The  theory  proposed  by  Perzyna  in  1963  [23-26]  generalizes  the  rigid/ perfectly-viscoplastic 
theory  of  Hohenemser-Prager  into  an  elastic/elastic-viscoplastic  theory  that  includes  hardening.  The 
theory  in  its  initial  form  is  for  small  deformation.  Prior  to  yield  the  deformation  is  considered  to  be 
rate-independent  and  recoverable,  i.e.,  elastic.  Whatever  inelastic  deformation  that  may  actually 
accompany  the  elastic  deformation  is  considered  to  be  insignificant. 


An  initial  ‘static’  yield  condition  is  assumed,  defined  by  the  yield  function 


FCT  E"  R 1  —  /(T’  E  )  . 

F(T,E  **)  *(>') 


in  which  the  possibility  for  isotropic  and  anisotropic  hardening  are  included  by  virtue  of  the  presence 
of  E"  and  R  as  variables.  The  existence  of  a  static  yield  condition  is  an  idealization  of  the  observa¬ 
tion  that  below  strain  rates  that  are  very  small,  a  strain  rate  effect  upon  yield  cannot  be  ordinarily 
observed. 


A  cornerstone  of  rate-independent  plastic  constitutive  theory  is  Drucker’s  postulate  of  positive 
plastic  work  [1,2],  which  leads  to  the  conclusions  that:  (i)  The  yield  surface  is  a  convex  surface  in 
stress  space,  and  (ii)  the  plastic  strain  increment  dZ",  or  [he  plastic  strain  rate  E",  must  be  normal 
to  the  yield  surface  F(T)  =0  at  all  points.  Thus  d E”  or  E"  must  also  be  parallel  to  the  gradient  of 
the  yield  function. 


dZ”  =A~, 
3T 


where  A  is  an  arbitrary  positive  scalar  function  having  the  same  arguments  as  the  yield  function. 


Considering  rate-independent  plasticity  as  the  lower  limit  of  a  rate-dependent  plasticity,  or  visco- 
plasticity,  the  inelastic  strain  tensor  is  postulated  to  be  normal  at  each  point  to  the  ‘static’  yield  sur¬ 
face,  and  to  all  subsequent  rate  dependent  yield  surface.  This  implies  the  relation 


E"  =  A- 


The  non-negative  scalar  valued  function  A,  other  than  possibly  being  a  function  of  the  same  variables 
as  F ,  is  arbitrary  otherwise.  Perzyna  takes  A  to  be  an  unspecified  functional  of  the  yield  function  F 
such  that 


A  =  7q  <<*>  (F)>  = 


$(F)  F>  0 
0  F<>  0 


where  70  is  a  constant.  It  is  presumed  that  the  functional  4>(F)  may  be  chosen  to  represent  the  results 
of  tests  on  the  behavior  of  metals  under  dynamic  loading.  With  the  choice  of  A  in  this  form 

_  dF  . ,  ^  3/  ,, 


E”  =  70  <W>  =  7 «W>  J£. 


The  strain  rate  is  assumed  to  be  decomposable  into  a  linear  fat  SIGMA  of  elastic  and  inelastic 
components 

E  =  E'  +  E'\  (4.6) 


K5«i 


(4.7) 


with  the  elastic  strain  rate  specified  by 

1  a  0  I  -  2v 

£'  =  —-<  D  +  ('rT)  1 

2/t  or  3£ 

The  rate  dependence  of  the  yield  function  is  not  obvious  from  Eq.  (4. 1).  However  from  (4.5) 

tr(±"  i ")  =  II' f.  =72*2  (F)fr(|£-  Jf) 


from  which 


*(F)  =  -VF“  j 

7 


•1/2 


E”  i  dT  dT 


Taking  the  inverse  of 


f  -  i  - 1  -  *“ 


V//V 

L 

K..K 

-1/2 

7 

L 

dT  dT 

V.  J 

J 

so  that 


fC T,E”)  =  £(**"') 


1 

^  H't"  f 

-  J  tr 

rl 

7  1 

dT  dT 

^  J 

j  j 

(4.8) 


The  dependence  of  the  yield  function  /  on  the  inelastic  strain  rate  appears  thru  the  dependence  of 
<J>_1  on  the  second  invariant  of  E”. 


In  summary: 

t  -  t'  +  t" 

E'  ,"T)1 
E”  -  7  <  *(f )  >  jjL  for  F  >0 
E”  «  0  for  F  SO 


M,  v,£,  70,  are  constants. 


(4.9) 


Comments: 

1.  Up  to  this  stage  of  development  of  the  theory,  the  assumption  of  inelastic  incompressibility  has 
not  been  introduced,  since  the  stress  tensor  rather  than  its  deviator  appears  as  an  argument  for 
the  yield  function.  Also,  the  specific  form  of  the  hardening  is  left  open. 

2.  As  with  the  Bingham  and  the  Hohenemser-Prager  theories,  inelastic  deformation  occurs  only  for 
F>0  rather  than  when  F  a 0,  as  in  the  rate-independent  theories  of  plasticity.  In  the  Bingham 
and  the  Hohenemser-Prager  constitutive  theories  the  functional  $(£)  have  the  specific  forms 


*(F)  =  f  =  l  - 


yJlFo  ' 

V  T 


respectively. 


3.  The  dependence  of  yield  upon  the  strain  rate,  as  indicated  by  Eq.  (4.8),  applies  for  subsequent 
yield  only.  The  initial  yield,  however,  in  this  theory  is  rate-independent.  To  see  this,  for  onset 

of  yield  F  =0,  and  as  <t(0)  =  0—  E"  =0  and  W"  =0.  Thus  -  1  =0,  so  that 

/( T,0)  -R{ 0)  =  constant,  which  represents  the  size  of  the  yield  surface  at  onset  of  yield, 
independent  of  the  strain  rate  at  which  the  yield  initiates.  For  uniaxial  stress  this  constant 
represents  the  yield  stress  Yt  or  Ys . 

4a.  Inelastic  Incompressibility  and  Isotropic  Hardening 

More  specialized  forms  of  the  constitutive  equations  can  be  written  as  additional  assumptions  are 
introduced.  If  inelastic  incompressibility  is  assumed  and  hardening  is  limited  to  isotropic  hardening 
only,  the  yield  function  (4. 1)  appropriate  to  these  conditions  has  the  form 


F(T,R) 


/(//'*,///'.) 

T  T 

R(W") 


(4.10) 


If,  furthermore,  for  the  function  /  a  von  Mises  form  is  adopted,  that  is,  /  =  (//'„)l/2 

T 


F(I,R)  ~  ' ' 


(4.11) 


and  Equations  (4.5)  and  (4.8)  accordingly  become 


E”  -  t  <* 


yflF. 


r  ’Tv^' 


(4.12) 


J1FI=R{W")  i+*-‘ 

v  T 


yfir~. 


(4.13) 


4b.  Uniaxial  Constitutive  Equations 

For  axial  loadings,  or  uniaxial  stress,  Eqs.  (4.9)  reduce  to  uniaxial  form  (cf.  Appendix  for  the 
stress  and  strain  tensor  reductions).  Assuming  inelastic  incompressibility,  isotropic  hardening  only 
and  a  von  Mises  form  for  the  function  /,  i.e.,  Eq.  (4.11)  for  the  yield  function,  a  reduction  to  uniax¬ 
ial  stress  gives,  in  summary  form. 


i  ■*!  >  tit  i.< 


r. 


€  =  €'+€" 


e  = T 


-Vf  7  <  ♦(F)  >  for  F  >  0 


«”  =  o  for  F  ■& 0 


(4.14) 


Vf 

W"  =  j  ode",E  and  7C  are  constants. 


Eq.  (4. 13)  correspondingly  reduces  to 


V?- 


/?(1P")  1  +  ♦ 


V?«' 


(4.15) 


The  following  forms  for  the  functional  ♦(F)  are  suggested  by  Perzyna  as  forms  that  may  possi¬ 
bly  be  used  to  correlate  uniaxial  stress-strain  experimental  data  for  structural  metals: 


♦(F)  =  F 
♦(F)  =  F" 

♦(F)  =  L  An  F" 

♦(F)  *  eF  -  1 

♦(F)  =  E  /4„  (er  -  1). 

n  —  1 


(4.16) 


5.  Chaboche  Theory 

The  viscoplastic  theory  proposed  by  Chaboche  [27-30]  is  fundamentally  similar  in  its  underlying 
structure  to  Perzyna’s  theory,  in  that  a  normality  assumption  of  the  infinitesimal  strain  rate  to  the 
yield  surface  defines  the  flow  rule.  In  this  theory  specific  forms  are  chosen  for  the  yield  function  F 
and  for  the  flow  rule  functional  $(F).  The  theory  also  introduces  specific  rate  equations  that  govern 
the  evolution  of  the  hardening,  which  is  taken  as  a  combination  of  isotropic  and  kinematic  types  of 
hardening. 

The  deformation  is  assumed  to  be  of  small  order  and  isothermal,  with  the  inelastic  part  con¬ 
sidered  as  being  volume  preserving.  The  yield  function  is  taken  to  have  the  von  Mises  form,  modi¬ 

fied  for  the  inclusion  of  hardening.  Thus 

F(T,Y,R)  =  F(T-Y)-R  -  Jll’o  ,  -R,  (5.1) 

v  cr-Y) 

where 

[O  O  0  0  1 

(T-Y)-(T-Y)J  ,  (5.2) 

(T-Y) 


Y  is  the  deviator  of  the  kinematic  hardening  tensor  and  R,  the  isotropic  hardening  variable,  is 
assumed  to  be  a  function  of  the  inelastic  strain  accumulation  p ,  defined  by  Eqs.  (1.10)  and  (ii.7)2. 


The  normality  assumption  [cf.  Eq.  (4.2)],  here  written  for  the  inelastic  strain  rate  rather  than  as 
its  increment,  requires,  in  view  of  Eqs.  (5.2)  and  (ii.  12)  that  in  the  space  of  the  stress  deviator 


*•  -  *4 

3T 


T-Y 

Jir .  . 

V  (T-Y) 


(5.3) 


The  scalar  proportionality  function  A  which  is  positive,  and  which  can  depend  possibly  on  the  same 
variables  as  F,  is  arbitrary  otherwise.  As  in  the  Perzyna  theory,  [cf.  Eq.  (4.4)],  it  is  chosen  as  a 
functional  of  the  yield  function  having  the  particular  power  form: 


A  -  <*(F)> 


♦(F)for  F>0 
0  for  F sO 


(5.4) 


where 


*(F)  *  (“■)" 


(5.5) 


and  K',n  are  temperature  dependent  material  parameters.  The  so-called  flow  rule,  comparable  to  the 
flow  rule  (4.5)  of  the  Perzyna  theory,  in  this  theory  has  the  explicit  form 


E” 


J1F77  -R 

y  (T-Y) 
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O  O 

(T-Y) 

Jir .  . 

y  (T-Y) 


for  F  >  0. 


I 

(5.6) 


The  hardening  that  accompanies  the  inelastic  deformation  is  described  by  a  pair  of  rate 


(5.7) 


equations.  For  the  isotropic  hardening  variable  R(P) 

-  B(Q -R)  ^  *  B(Q  - R)  |  t"  ! 
dt  at 


where 

P  -  f  d  |  E"  |  =  P 

to  Jo  dr  L 

dr 

(5.8) 

since 

d  |E"|  [rr^(E”-E">]  = 

d(yJlI'r.). 

(5.9) 

It  follows  therefore  that 

P  *  \H"  | 

(5.10) 

For  the  anisotropic  part  of  the  hardening,  the  kinematic  hardening  that  is  employed  is  described  by 
the  rate  equation 


dY 

dt 
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-  r (//',) 
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(5.11) 


In  these  hardening  rate  equations,  B,  Q,  C,  A,  I\  and  m  are  material  parameters  that  are  temperature 
(frpfpA»nt  The  first  contribution  to  the  kinematic  hardening  rate  equation  is  Prager’s  form  [5]  which 
is  appropriate  for  monotonic  loading.  The  second  term  generalizes  this  hardening  rule  so  as  to  allow 
it  to  apply  for  change  of  load  direction,  or  for  cyclic  loading  [32].  The  third  models  the  hardening 
recovery,  or  softening,  that  is  often  observed  at  high  temperatures.  The  isotropic  hardening  rate 
first-order  differential  equation  (5.7)  can  be  solved  directly,  giving 

R(P)  -  Q+iRo-Q*-” .  (5.12) 


As  the  iM-rnmnlatrH  inelastic  strain  P  —  0,  R  (0)  —  R*,  which  represents  the  initial  yield  stress.  As 
the  iMfimmlatwt  inelastic  strain  becomes  large,  R  —  Q,  which  represents  the  ‘saturation’  value  for  R. 
In  the  absence  of  elevated  temperature  the  hardening  recovery  (or  thermal  softening)  term  in  Eq. 
(5.11)  is  dropped,  i.e.,  T  -  0. 


In  summary: 

E  -  E'  +  E”. 
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dR 

dt 

R 

O 

dY 

dt 
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0  for  FsO. 
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(T-Y) 
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(T-Y) 


JTF7 T  -  R(P). 

V  (T-Y) 

B(Q-R .  p  -  |E”|. 


-  Q+(R0 -Q)e-". 

»  C(AH"  -  VjY  I  E”  !)-r(//„) 


for  F> 0. 


Y. 


(5.13) 


\n  ,v,E\  are  elastic  constants  ,  \K' ,n  ,B,  Q,A  ,  C,r,m,Ra  )  are  viscoplastic  material  parameters. 


5a.  Uniaxial  Constitutive  Equations 


Referring  to  the  Appendix  for  the  reductions  and  notational  changes  adopted  for  axial  loading, 
Equations  (5.13)  may  thereby  be  reduced  to  uniaxial  form.  In  so  doing,  the  following  transformations 
are  also  utilized 


yfip-p  |  de"  |  =  j'|e"|  dr,  p 

yf^R(P)-R(p)  ,  .  (j)~r^7 

u  + 1  *  ■ 

<{)  2  v  JA~a  .  7,,-* 

\[jQ-q,  VjC-c,  yu-y. 


where  «"  =  £[,”  and 


For  unaxial  loads  the  constitutive  equations  corresponding  to  Eqs.  (5.13)  can  be  expressed  as  follows: 
e  *e'  +  is" 
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ft 
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*0  forFjsO. 


Sg{a-Y) 
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for  F  >  0. 


f  -  |n-y|  -R(p) 

=b{q-R)^  ,  p=  |i"|  (5.14) 

R  m  q  +(R0-q)e~bp ' 

~  “C(a«"  -  K  |«"|)-tt|  K|mSg(lr)  ,  Sg(K)  - =  ±1 
(Ar,n,6,q,c,a,7,m,R0}  are  viscoplastic  material  parameters. 


The  non-linear  first-order  differential  equation  for  the  kinematic  hardening,  Eq.  (5.14)g,  must  be 
solved  numerically.  In  circumstances  for  which  the  hardening  recovery  term  can  be  neglected,  i.e. 
7*0,  the  remaining  differential  equation,  which  is  homogeneous  in  time,  can  be  brought  to  a 
separated  form  and  integrated  directly,  giving 

r  -  so  +(y0-ttj)e"'c“"''"’),  (5.15) 


where  s  *  +  1  for  de”> 0,  s  *  -1  for  de"  <  0,  and  Y0  and  e"0  represent  the  fixed  values  of  Y 

and  «"  at  the  previous  reversal  of  direction  of  the  strain  (or  load)  for  non-monotonic  and  cyclic  load¬ 
ing.  For  strictly  monotonic  loading 

Y  =a(l  -e~ct").  (5.16) 


Inversion  of  Eq.  (S.14)  for  the  axial  stress  yields 


a  =  Y  +  R  +K(e")n  . 


(5.17) 


For  monotonic  loading  from  zero  stress 


a  =  a(l  -  e~ct  )  +  q  +(R0-q)e 


+  *(e'')"  • 


(5.18) 


5b.  Modification  for  the  Strain  Rate  Dependence  of  Initial  Yield 

As  in  Perzyna’s  theory,  the  Chaboche  theory  also  fails  to  account  for  the  fact  that  the  onset  of 
yield,  or  the  initial  yield,  is  strain  rate  dependent.  When  yield  commences  e"  =  e"  —  p  =0,  and 
from  Eq.  (5.14),  lim  R(p)  -R0,  whereas  from  Eq.  (5.18)  lim  a  =R0.  Thus  R0  represents  the 

P-  0  («",i'')-0 

initial  yield  stress  for  either  uniaxial  tension  or  compression,  the  value  of  which  changes  with  the  rate 
of  strain  e  at  which  it  takes  place.  To  incorporate  this  dependency  into  the  constitutive  equation,  the 
isotropic  hardening  variable  can  be  taken  to  be  a  function  of  both  the  accumulated  inelastic  strain  p 
and  the  strain  rate  e,  where  the  dependence  on  e  comes  through  its  appearance  as  an  argument  for  the 
initial  yield  stress  R0.  The  Chaboche  theory  can  be  modified  accordingly  by  assuming  that  the  isotro¬ 
pic  hardening  variable  has  the  form  [32-34] 

R(p ,  e)  =  q  +  [/?<>(«>  ~  qV  (5- 19) 

where 

lim  R(p,k)  -  Rq(0)  »  Roo  (5.20) 

o 


is  the  initial  yield  stress  at  quasi-static  strain  rate,  that  is,  the  strain  rate  below  which  no  significant 
observable  effect  upon  the  yield  stress  can  be  detected  by  the  usual  means  of  measurement. 


The  form  of  the  initial  yield  function  Rrfe)  can  be  determined  by  or  suggested  from  initial  yield 
stress  versus  strain  rate  test  data  covering  the  strain  rate  range  of  interest.  As  an  illustration,  test  data 
gathered  for  the  alloy  Inconel  at  1200°F  indicates  that  the  dependence  of  the  initial  yield  stress  on  the 
strain  rate  is  significant  only  for  the  low  strain  rate  range  <  <  10~3  sec  -l  [35-38].  From  the  experi¬ 
mental  information  provided  by  these  sources  it  was  possible  to  correlate  the  initial  yield  dependence 
upon  strain  rate  over  the  range  10~9<e<  10-3  sec  -1  by  means  of  the  relation  [32,33,63] 


R  o(f)  *  Roo 


1  +  ln(l  +  exp 
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E  ctj  (In  e)' 
i  «o 


(5.21) 


In  this  expression  the  a,  are  the  correlation  coefficients  that  are  obtained  from  non-linear  regression 
aqalysis  of  [he  data.  This  form  can  be  generalized  to  multi-dimensional  stress  by  replacing  e  by 
I  El  -[fr(E  E)]‘/2. 


PART  n 

VISCOPLASTIC  THEORIES  WITHOUT  A  YIELD  CONDITION 


6.  Bodner-Partom  Theory 

The  initial  formulation  of  the  Bodner-Partom  theory  for  viscoplastic  behavior  [39-41]  incor¬ 
porates  only  the  isotropic  contribution  to  the  inelastic  hardening.  The  theory  is  constructed  within  the 
framework  of  small  isothermal  deformation  with  the  infinitesimal  strain  divided  into  elastic  and  inelas¬ 
tic  parts.  Thus 

E  =  E'  +  E",  (6.1) 

with  E'  the  linear  elastic  strain  rate  given  by  Eq.  (4.7).  The  inelastic  component  of  the  deformation 

O 

is  assumed  to  be  incompressible  so  that  E"  =  E".  With  the  absence  of  a  yield  condition,  inelastic 
strain  always  takes  place  at  all  levels  of  stress,  however  small. 

The  deviator  of  the  inelastic  strain  rate  is  assumed  to  be  proportional  to  the  stress  deviator 

4~  (E”)  =~(E'')  -  E”  **T  (6.2) 

at  at 

where 

*»*<//'.,Z)  (6.3) 

T 

and  Z  is  an  unspecified  scalar  internal  variable  that  is  determined  by  a  rate  equation  having  the  form 

^  »/(T,E”,Z).  (6.3a) 

at 

It  follows  from  Eq.  (6.2)  that 

tr(E''-  E”)  -  Hi...  *  *V(T  -  T)  *  // '»*2  (6.3b) 
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from  which 


and 


E”  -VFZ 


o 
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V"V 


(6.4) 


(6.5) 


The  development  leading  from  Eqs.  (6.1)  to  (6.5)  essentially  parallels  the  development  of  the 
Saint  Venant-Levy  theory  for  perfectly  plastic  flow,  except  that  the  scalar  internal  variable  Z  has  been 
introduced  to  serve  as  the  isotropic  hardening  variable.  To  obtain  a  specific  relation  between  the  har¬ 
dening  variable  and  the  inelastic  strain  rate,  recourse  is  made  to  the  Gilman-Johnston  microstructural 
model  for  the  motion  of  dislocations  in  crystalline  solids.  According  to  this  model  the  speed,  V ,  of  a 
moving  dislocation  is  related  to  the  resolved  shear  stress,  r,  by  the  relation 


In  this  equation  v*  represents  the  limit  value  of  v  as  r—  oo ,  and  thus  represents  the  limit  value  of  the 
dislocation  velocity.  The  term  d  is  called  the  ‘drag  stress’  that  impedes  the  motion  of  the  dislocation, 
i.e.,  in  macroscopic  terms,  the  source  of  the  hardening.  Bodner-Partom  adopt  this  model,  generaliz¬ 
ing  it  to  multidimensional  stress  and  setting  it  to  a  macroscopic  description  by  means  of  the  following 
correspondences: 


T-U'o  ,  d~Z  ,  V*—Dq  , 
T 


and  the  postulated  form 


Dq  exp  - 


The  internal  variable  Z  is  thus  identified  with  the  ‘drag  stress’  as  a  measure  of  the  isotropic  hardening 
associated  with  the  inelastic  deformation.  The  parameter  D0  is  interpreted  as  the  limit  value  of  the 
shear  strain  rate.  Combining  Eqs.  (6.5)  and  (6.7)  gives 


E”  -  D0  exp 


as  the  flow  rule,  which  includes  Z  explicitly  as  the  isotropic  hardening  variable.  The  constants  D0 
and  n  are  treated  as  material  parameters  that  are  temperature  dependent. 

The  evolution  of  Z,  that  is,  its  growth  with  inelastic  deformation,  is  described  by  the  postulated 
rate  equation 

?-<*■-*  >¥- 


which  has  the  solution 


zi  +  (z0-zx)e' 


(6.10) 


The  quantities  m ,  z0  and  z ,  are  material  parameters,  where  z0  is  the  value  of  Z  when  the  inelastic 
strain  is  zero.  Note  the  similarity  in  the  form  of  the  isotropic  hardening  Eq.  (6. 10)  to  that  of  Eq. 
(5.12)  of  the  Chaboche  theory.  In  Eq.  (6.10)  the  hardening  is  expressed  in  terms  of  the  inelastic 
work,  whereas  in  Eq.  (5. 12)  it  is  determined  by  the  accumulated  inelastic  strain.  The  parameter  z , 
represents  the  ‘saturation’  value  of  Z  as  W  —  oo,  and  corresponds  in  this  sense  to  the  parameter  Q 
in  Eq.  (5.12).  Also  note  that  while  Rq  in  Eq.  (5.12)  can  be  identified  as  the  initial  yield  stress,  no 
such  interpretation  can  be  given  to  the  parameter  Z0.  In  circumstances  where  hardening  recovery 


takes  place  (e.g.  thermal  softening  at  elevated  temperature)  the  isotropic  hardening  rate  equation  (6.9) 
is  modified  by  the  addition  of  hardening  recovery  term. 


dZ  .  dW" 

-j  --»(z,-Z)  —  -At, 


in  which  Atz2  and  r  are  additional  material  parameters. 

In  summary: 

E  »  t'+  E"- 

•  1  3  0  1  —  2i» 


•7  V 

2  ~z2 
Zi 


(6.11) 


E"  *  D0  exp  - 


dZ  .  dW 
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(6.12) 


2  *  zi  +(z0“Zi)«  m,r 

(m,  »>,£}  are  elastic  constants  ,  (D0,  n ,  z0,  z  j ,  z2,  m ,  A ,  r )  are  viscoplastic  material  parameters. 

6a.  Modification  for  Directional  Hardening:  Stouffer-Bodner  Theory 

A  modification  of  the  Bodner-Partom  theory  was  proposed  by  Stouffer  and  Bodner  [42,43]  that 
allows  it  to  incorporate  the  Bauschinger  effect  and  the  generally  anisotropic  nature  of  plastic  harden¬ 
ing  associated  with  load  reversals.  In  the  absence  of  a  yield  surface  to  provide  a  geometrical 
representation  of  hardening,  the  anisotropic  hardening  that  is  introduced  is  referred  to  as  a  ‘directional 
hardening’  that  appears  to  be  similar  to  the  more  familiar  kinematic  hardening. 


Equation  (6.2)  is  replaced  by  the  relation 


£(=”>-£(  e")-a-t. 


(6.13) 


in  which  A  is  a  fourth-order  tensor  which  governs  the  anisotropy  of  the  inelastic  behavior.  This 

O 

equation  implies  inelastic  incompressibility  since  E"  3  E.  Because  of  the  symmetry  of  the  tensors 

.  o 

E"  and  T,  A  must  have  the  symmetries 

*tjki  *  A jju  *  A ,  (6.14) 

which  reduces  the  number  of  independent  non-zero  components  of  the  tensor  from  81  to  36.  This 
symmetry  allows  introduction  of  the  following  (rotational  change,  where  the  components  of  the  tensors 

T  and  E"  are  re-labeled  according  to  the  scheme 


Th— T|  , 

2T|2— T4 

O  0 

T22-T2  , 

O  O 

2T23—T5. 

(6.15) 
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(6.16) 


1  m.'M  m't  .fa  +_>  |,|  |  »>l  1-, 


ll.  •< 


Expressions  (6. 14)  and  (6.  IS)  enable  Eq.  (6. 13)  to  be  expressed  in  the  form 

6  o 

sa"=  E  A^,  a  =1,2 . 6.  (6.16) 

(3-1 

The  independent  non-zero  components  of  the  fourth-order  tensor  A  are  thus  represented  by  the  6  x  6 
matrix  [A^.  This  matrix  is  next  assumed  to  be  symmetric,  although  no  explanation  is  given  as  to 
what  this  means  or  implies  in  terms  of  reduction  of  the  degree  of  anisotropy.  (Analytically  speaking 
the  assumed  symmetry  reduces  the  number  of  independent  components  from  thirty-six  to  twenty-one). 
With  the  symmetry,  the  principle  values  (eigenvalues),  A^,  a  =  1,. . .  ,6  of  [A^J  are  real  and  the 
principal  directions  (eigenvectors),  ea,  are  orthogonal.  Knowledge  of  the  eigenvectors  permits  con¬ 
struction  of  an  orthogonal  matrix  (tensor),  Q~‘  =  Qr,  such  that  by  means  of  the  similarity  transfor¬ 
mation  Q-  A  •  Qr,  reduces  A  to  its  diagonal  form 


[Aa05<Xt]  ,ej 


(6.17) 


Relative  to  the  principal  direction  basis  ea,  a  =  1 , ....  6,  the  constitutive  equations  (6. 16)  take  the 
form 


£a”  -  E  Agg  Tff 
0*1 


(6.18) 


Aaa>Ct  =  0 


a  =  1 , . . . ,  6. 


0,  ant/8 


From  Eq.  (6. 17)  it  follows  that 


/ c"  *  E"i  +  £”2  +  £"3  =  AjjT]  +  A22T2  +  A33T3. 


(6.18) 


For  isotropic  hardening,  the  hardening  tensor  can  have  only  two  independent  non-zero  components. 
Thus  A  n  »  A  22  m  A  33  *>  X  can  be  one  of  them  and  ,  consequently 


OOO 


/j>i  *  X  (Tj  +  T2  +  T3)  —  \  la  =0 


for  all  values  of  X.  This  implies  plastic  incompressibility.  For  anisotropic  hardening  in  general. 
An  *  A 22*  A  33  and  /j.,  *0,  which  implies  plastic  compressibility  in  contradiction  with  the  inelastic 
incompressibility  implied  in  Eq.  (6.13). 

Paralleling  the  assumption  (6.3)  for  isotropic  hardening,  the  principle  values  of  the  anisotropic 
material  response  tensor  A  are  assumed  to  be  functions  of  the  second  invariant  of  the  strain  deviator 
tensor  and  of  a  tensorial  internal  variable  Z,  where  the  components  of  Z,  referred  to  principal 
directions  of  the  matrix  [A^],  are  such  that  »  0,  a  *  0.  Thus 


Fa(ir  o,Zm) ,  a  =  1,...,6  (no  sum  on  a). 
T 


(6.19) 


The  internal  variable  Z  is  interpreted  as  a  hardness  variable  that  represents  the  resistance  to  anisotro¬ 
pic  inelastic  deformation.  The  form  of  Eqs.  (6.19)  is  taken  to  be  similar  to  the  form  of  Eq.  (6.8). 
That  is,  Stouffer  and  Bodner  ascribe  to  anisotropic  hardening  the  same  microstructural  ‘drag  stress' 
mechanism  modeled  by  Eq.  (6.6),  that  they  used  to  introduce  isotropic  hardening  into  the  Bodner- 
Partom  theory  when  they  assume  that 


Aa „  =  D0  exp 
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a  =  1 , ....  6  (  no  sum  on  a) .  (6.20) 


The  hardness  internal  variable  Z  is  assumed  to  be  governed  by  a  rate  equation  similar  to  the  rate 
equation  (6.11)  for  isotropic  hardening.  To  introduce  directionality  and  the  Bauschinger  effect,  a 
material  parameter  q  is  introduced  that  serves  to  describe  the  relative  amount  of  the  isotropic  part  to 
the  directional  part  of  the  hardening  in  the  following  manner: 

Zw  =  qZ+(l-q)Zua  for£”a>0 

a  =  1,. . .  ,6,(  no  sum  on  a)  ,  (6.21) 

Z*,  -  *Z-(1  -q)ZUa  for  L"a  >  0 


where 


Un 


S", 
i  in 


are  the  ‘direction  cosines’  of  the  rate  of  strain 


(6.22) 


and 


Z  *  m(z,  -Z)W"  -  Az | 


r_  ^ 
Z  -z2 


(6.23) 


Integration  of  Eqs.  (6.21)  gives  for  each  a  »  1,. . .  ,6  , 


Z*n  +  Q  jQZ(r)dr  +  (1  -q)Sg(£a")\Z(T)Ua(T)dT, 


Sg  (Ea") 
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l=."  I 


=  ±1. 


(6.24) 


A  value  of  q  ■  1  indicates  completely  isotropic  hardening,  whereas  q  -0  signifies  that  the  hardening  is 
entirely  directional. 


For  eqch  value  of  a 
strain  rate  Ea"  as 
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1,. . .  ,6,  expressions  (6.18)  and  (6.20)  give  the  six  components  for  the 
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in  which  the  six  directional  hardness  variables  Zm  are  obtained  by  means  of  Eqs.  (6.22)  and  (6.23). 


In  summary: 
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(/*,  i»,  £ }  are  elastic  constants,  {Z>0»  « ,  Zo»  *  l >  * 2» m  •  ^ » r  •  ? )  ***  viscoplastic  material  parameters. 


Comments: 

1.  The  anisotropic  hardening  that  is  described  by  the  StoufFer-Bodner  theory  is  an  attempt  to  simu¬ 
late  the  kinematic  hardening  appearing  in  yield-based  theories  without  the  benefit  of  the 
geometrical  description  of  hardening  afforded  by  a  yield  surface.  The  anisotropic  hardening, 
described  by  a  9  x  9  matrix  in  the  most  general  sense,  is  reduced  by  symmetry  considerations 
and  assumptions  so  as  to  allow  the  main  diagonal  components  of  the  6  x  6  reduced  hardening 
matrix  to  be  taken  as  functions  of  the  principal  values  of  a  hardening  tensor  Z,  which  is  con¬ 
sidered  to  be  an  internal  variable.  The  components  of  the  hardening  tensor  are  further  assumed 
to  be  parallel  to  the  components  of  the  inelastic  strain  rate  tensor.  In  this  way  a  directional  har¬ 
dening  analog  to  kinematic  hardening  is  constructed. 

2.  The  functional  form  used  to  relate  the  components  of  the  material  tensor  A  to  the  com¬ 
ponents  Zaa  of  the  internal  variable  hardening  tensor  Z  is  identical  to  the  one  used  for  the  iso¬ 
tropic  hardening  theory.  In  other  words,  the  same  microstructural  model  is  used  to  suggest 
macroscopic  relations  that  describe  two  different  kinds  of  hardening,  although  isotropic  and 
anisotropic  hardening  are  generally  believed  to  be  caused  by  different  microstructural  mecha¬ 
nisms. 

3.  The  parameter  q  which  proportions  the  hardening  into  isotropic  and  directional  parts  is  arbi¬ 
trarily  assigned  a  value.  It  is  not  measured.  For  cyclic  loading,  where  during  the  course  of  the 
cycling  the  character  of  the  hardening  changes,  the  value  chosen  for  q  would  have  to  be  changed 
repeatedly  as  the  cycling  progresses. 

4.  The  parameter  D0,  representing  the  limit  value  of  the  shear  strain  rate,  is  also  a  parameter  that 
is  not  measured.  A  value  of  104  sec-1  was  initially  suggested  by  Bodner.  Subsequently  values 
of  10^  and  108  have  been  assigned  for  the  value  of  D0  by  others  using  this  viscoplastic  constitu¬ 
tive  theory. 


[D0,n,Zo,Z\,Zi,m,A,r ,q\  are  viscoplastic  material  parameters 


7.  Walker  Theory 


The  viscoplastic  constitutive  theory  proposed  by  Walker  that  is  shown  here  is  the  theory  as  of 
1981  [44,45].  Walker  advances  two  theories,  an  integral  version  and  a  rate  version.  Both  theories 
are  appropriate  for  small  isothermal  deformations  that  are  assumed  to  be  decomposable  into  elastic 
and  inelastic  parts. 

7a.  Integral  Theory 

The  integral  version  is  obtained  by  modification  of  the  standard  3-parameter  spring-dashpot 
model  for  a  viscoelastic  solid.  For  a  uniaxial  stress  it  may  be  shown  to  give  the  rate  equation 

a  +  X  a  =  Eft  +  X  E,t  ,  (7.1) 

where  Ef,  Es  and  X  are  the  three  material  parameters  of  the  model. 

The  parameter  Ef ,  representing  the  elastic  modulus  of  the  spring  in  series  with  the  elastic  spring  and 
dashpot  that  are  arranged  in  parallel  to  each  other,  gives  the  instantaneous  response  of  the  system 

a  —  Eft  (7.2) 

which  is  the  first  term  of  the  right  side  of  Eq.  (7.1).  For  very  slow  rates  of  loading  or  strain,  the 
response  of  the  model  is  given  essentially  by  the  second  term 

E,t  *  Y,  (7.3) 

which  is  referred  to  as  the  ‘equilibrium  stress’.  Es  is  the  elastic  modulus  of  the  spring  in  parallel  to 
the  dashpot.  The  parameter  X,  representing  the  action  of  the  dashpot,  is  interpreted  as  a  viscosity 
coefficient. 


a 


The  Boltzman  integral  form  of  the  standard  3-parameter  model  for  a  viscoelastic  solid  subject  of 
uniaxial  stress  is  given  by 


dr 


(7.4) 


where 


G(t)  =  Eft-*  +  £j(l  -  e~u) 


(7.5) 


is  the  relaxation  function.  Substitution  of  expressions  (7.5)  and  (7.3)  into  (7.4)  leaves 
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it  follows  that 
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(7.7) 


The  modification  of  the  above  viscoelastic  constitutive  theory  to  one  representing  viscoplasticity 
is  proposed  in  the  following  manner.  Inelastic  strain  is  defined  as  the  difference  between  the  total 
strain  and  the  elastic  strain 

t"  =  e  — jr o  ,  where  E  —  Ef .  (7.8) 

£ 

The  linear  elastic  ‘equilibrium  stress’  relation  (7.3)  is  replaced  by  a  nonlinear  relation  that  depends  on 
the  inelastic  strain  and  its  accumulation 

Y=Y(e”,p).  (7.9) 

The  viscosity  coefficient  X  is  assumed  to  be  nonlinearly  proportional  to  the  inelastic  strain  accumula¬ 
tion  rate 

F  0 — -) 

X  =  ~p  ",  (7.10) 


with  R  and  n  as  constants. 

In  view  of  Eq.  (7.10), 
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where 
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Substitution  of  expressions  (7.9)  and  (7.11)  into  Eq.  (7.7)  yields 


ait)  =  Y{t)  +  |o 
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Expression  (7.9)  for  Y  is  postulated  to  have  the  form 


Y{t)  =  Yq  +n  i  e"(t)  +  ~lZ(,)_Z(T)l  \  ~-\  dr 


in  which  Z(t)  is  specified  by 

Z(0  =  [«3  +  ntpit)]  In 
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dr 


(7.14) 


(7.15) 


Equations  (7.8)  and  (7.12)  -  (7.15)  constitute  Walker’s  proposed  viscoplastic  constitutive  theory. 


The  ‘equilibrium  stress’  Y  is  interpreted  as  an  internal  variable  representing  the  microstructural 
‘back  stress’  which  —  “is  a  measure  of  the  globally  averaged  dislocation  arrangement  and  produces 
directional  hardening’’.  The  nte''  contribution  to  Eq.  (7.14)  is  meant  to  model  the  Bauschinger 
effect.  Expression  (7.  IS)  is  intended  to  model  cyclic  hardening  and  softening,  where  the  logarithmic 
term  allows  for  cyclic  hardening  and  the  integral  term  for  the  decrease  of  the  equilibrium  stress  due  to 
thermal  recovery  (softening)  at  elevated  temperatures.  The  constant  R  in  Eq.  (7.12)  is  identified  as 
the  microstructural  ‘drag  stress’  usually  associated  with  isotropic  hardening,  [cf.  Eq.  (6.7)  of  the 
Bodner-Partom  theory].  It  is  here  taken  to  be  constant. 


In  summary: 
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Y(t)  =  YQ  +  nie"it)  +  (Qn2e-lZ(li-z(^\^-\dT 

Z(r)  »  [n3  +  n4/»]  In  ~4~p~  +  1  +  I0'*7 1  K  ~  d-i 


\n,R,Y0,m,n ,,. . .  ,n7]  are  temperature  dependent  viscoplastic  material  parameters. 


7b.  Rate  Theory 


A  viscoplastic  constitutive  theory  having  a  rate  form  is  also  constructed  from  the  3-parameter 
spring-dashpot  model  for  viscoelastic  material  response.  After  combining  Eqs.  (7.1)  and  (7.3)  and 
relabeling  Ef  as  E,  it  follows  that 


e  *  -p[ff  +  X(a-T)]. 

E 


(7.17) 


The  conversion  to  a  theory  for  viscoplasticity  begins  by  differentiating  Eq.  (7.8)  for  the  definition  of 
the  inelastic  strain 


(7.18) 


This  is  substituted  into  Eq.  (7.17),  along  with  Eq.  (7.10)  which  replaces  the  spring-dashpot  viscosity 
parameter  X  by  the  proposal  power  relation  of  the  inelastic  strain  rate.  The  result  is 
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Since  p  *  |  e"  | ,  Eq.  (7. 19)  may  also  be  expressed  by  the  equivalent  form 
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In  the  integral  theory  the  ‘equilibrium’  stress  Y  is  interpreted  as  an  internal  variable  that 
represents  the  microstructural  ‘back  stress’  that  is  associated  with  directional  hardening.  It  is  given 
by  the  proposed  relations  (7.14)  and  (7.15).  For  the  rate  theory,  Y  is  now  governed  by  the  proposed 
rate  or  evolution  equation 


where 


-J-  =(/»(+  n2)  t"  —  (Y  —  n  if")Z  , 


—  =  (n2  +  nie~n,p)p  +n6  |  Y|m_1 


(7.21) 


(7.22) 


In  the  integral  theory  the  constant  R  appearing  in  Eq.  (7.12)  was  interpreted  as  representing  the 
microstructural  ‘drag  stress’.  Here  it  is  assumed  to  be  a  variable  (an  internal  variable)  that  accounts 
for  isotropic  hardening  according  to  the  proposed  relation 


R  =  Ri-R2e~niP 


(7.23) 


In  summary 


e  *€-«=€  — -a 
E 


Sg(a-Y) 


R  ~Rl-R2e~'hP 
“  -  (n,  +n2)'e"  —  (Y  —nl€,')Z 

»(n3  +  n4e-">;,)p+n6|r|— 1 


(7.24) 


{n ,m ,nl,n2,n2,n4,n5,n6,n1,Rl,R2\  are  temperature  dependent  viscoplastic  material  parameters 


Comments: 


Generalization  of  Eqs.  (7.16)  and  (7.24)  to  multiaxial  states  of  stress  is  straightforward. 

Although  Walker  claims  that  the  two  theories  are  equivalent,  i.e.,  that  the  rate  theory  may  be 
obtained  from  the  integral  theory  by  differentiation  of  the  latter,  it  is  difficult  to  see  how  this  is 
in  fact  possible. 

The  development  of  both  of  these  theories  have  the  earmarks  of  an  ad-hoc  piecemeal  construc¬ 
tion. 

The  similarity  of  several  of  Walker’s  rate  equations  with  those  of  the  Chaboche  theory  should  be 
noted.  In  particular  Eq.  (7.24)  with  Eq.  (5. 14)3  for  the  inelastic  strain  rate;  Eq.  (7.24)3  with 
Eq.  (5.14)7  for  the  modeling  of  the  isotropic  hardening;  and  Eqs.  (7.24)45  with  Eq.  (5.14)g  for 
the  directional,  i.e.,  kinematic  hardening. 


8.  Kremple-Liu  Theory 

The  Kremple-Liu  theory  [46,47]  although  lacking  a  yield  condition,  is  nevertheless  constructed 
such  that  for  loading,  or  unloading  from  a  plastic  state,  the  predicted  model  response  will  be  elastic. 
This  theory,  for  small  isothermal  deformation,  is  also  developed  by  modification  of  a  3-parameter 
spring-dashpot  model  for  a  solid  uniaxial ly  loaded.  Starting  with  the  rate  equation  derived  from  the 
spring-dashpot  model, 

-£-(£,  +  £„)«  +  £.«  =  +a  (8.1) 

£j  E, 

in  which  £, ,  £„  are  the  elastic  coefficients  and  X  the  viscous  coefficient  of  the  model,  a  viscoplastic 
theory  is  proposed  by  changing  the  material  coefficients  from  constants  to  functions  of  the  stress  and 
strain. 

An  equilibrium  (zero  strain  rate)  stress-strain  relation  is  assumed 

a  =  4(e) .  (8.2) 


An  overstress  a  —  4(e)  is  defined  (similar  to  what  appears  in  the  Perzyna  and  Chaboche  theories). 
The  material  parameters  appearing  in  the  rate  equation  (8.1)  are  replaced  by  unspecified  material 
functional  of  the  overstress: 

~~(ES  +Em)~*{o  -4(e)) 

X(a-4(e)), 

while  the  stress  £,.e  is  replaced  by  the  equilibrium  stress  4(e).  The  proposed  viscoplastic  constitu¬ 
tive  equation  for  uniaxial  stress  is  thus 

♦(a  -  4(e))e  +  4(e)  =  X{o  -  4(e))ff  +  a ,  (8.3) 

in  which  the  material  functional  4  and  X,  as  well  as  the  function  4,  are  to  be  determined  by  experi¬ 
ment.  These  functions  are  required  to  be  bounded  continuous  functions  that  are  also  positive  and 
even.  They  are  also  required  to  be  such  that  for  all  values  of  their  argument  the  ratio  of  4  over  X  is 
constant,  equal  to  the  value  of  the  elastic  modulus. 

3* (8.4) 

X(<7  -  4(e)) 


The  form  of  Eq.  (8.3)  is  such  that  for  quasi-static  stress  or  strain  rates,  e  and  a  are  negligibly  small 
leaving  the  equilibrium  stress-strain  relation  (8.2),  while  for  very  high  stress  or  strain  rates  condition 
(8.4)  ensures  a  linear  elastic  response. 

Equations  (8.3)  and  (8.4)  together  may  be  shown  to  produce  the  following  additional  qualitative 
features  of  the  proposed  theory. 

(i)  There  is  initial  linear  elastic  behavior  for  all  rates  of  loading. 

(ii)  An  elastic  slope  is  obtained  upon  large  instantaneous  changes  in  strain  (or  stress)  rates,  e.g.,  at 
load  reversals  from  the  plastic  range. 


(iii)  Stress-strain  curves  obtained  at  different  rates  of  loading  are  non-linearly  spaced. 


(iv)  In  predicting  a  creep  test  with  tensile  creep  stress,  the  model  yields  e  >  Q  and  e  <  0  for  the 
total  strain  rate,  provided  that  (d4>/de)  a  0.  Primary  and  secondary  creep  can  be  reproduced. 


(v)  In  predicting  a  stress  relaxation  test  for  positive  hold  strain,  the  model  calculates  a  <  0  and 
a  2:  0,  where  relaxation  terminates  at  the  equilibrium  stress  a  =  <f>(e). 


Equations  (8.3)  and  (8.4),  when  combined,  give 

<--4  + . , 

E  EX( r  -  $(£)) 


(8.5) 


Since  the  first  term  is  the  elastic  stress  rate,  the  second  term  is  interpreted  as  the  inelastic  strain  rate. 
Thus 

«*£'+£"  (8.6) 


where 


a  -  $(<) 
EX  (a  ~  *(«)) 


(8.7) 


is  defined  as  the  viscoplastic  flow  rule  of  the  theory.  The  equilibrium  stress  function  <f>  and  the 
material  response  functional  X  for  a  given  material,  are  determined  from  uniaxial  test  data  at  the  tem¬ 
perature  of  interest. 


An  example  of  particular  representations  for  these  functions,  taken  from  references  [48,49],  are 
as  follows: 


$(«) 


E,(t-t o) 


X 

+ 


E-E, 

—  log 

2Ctanh(CX-3)  8 


cosh  U 
cosh  V 


+  ®o. 


(8.8) 

(8.9) 


where 

U  -  C[X  +  (£  -  to)]  -  3 

V  -  C[*  -  (e  -  eo)]  -  3  (8.10) 

C  *  3.63/X . 


E,  is  a  tangent  modulus  and  A,B,C,  are  parameters  determined  from  uniaxial  tests.  The  stress  and 
strain  values  a0  and  (q  are  values  at  points  where  a  load  reversal  occurs. 

Comments: 

1 .  The  above  theory  has  been  used  to  model  the  room  temperature  uniaxial  stress  behavior  of  AISI 
304  stainless  steel,  featuring  changes  in  strain  and  stress,  instantaneous  changes  in  strain  and 
stress  rates,  creep  and  stress  relaxation. 

2.  A  generalization  to  a  multi-dimensional  form  of  Eqs.  (8.3)  -  (8.5)  will  give: 

*(T  -  *(E))E  +  +(E)  -  X(T  -  *(E))T  +T .  (8.11) 

(8.12) 
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and 


9.  Miller  Theory 


This  constitutive  theory  is  developed  for  small  isothermal  deformation  and  uniaxial  stress.  It  is 
constructed  to  model  monotonic  and  cyclic  hardening  and  creep  behavior.  The  inelastic  component  of 
the  total  strain  rate  is  assumed  initially  to  have  the  form 


e 


»  t 


=  * 


a-Y 

R 


(9.1) 


where  ♦  is  an  unspecified  functional  of  the  argument  shown.  The  variables  Y  and  R  are  referred  to, 
respectively,  as  the  ‘rest  stress’  associated  with  kinematic  hardening,  and  the  ‘drag  stress’  associated 
with  isotropic  hardening.  To  give  specific  form  to  the  functional  *,  Miller  considers  the  following 
empirical  relation  which  relates  steady  state  creep  strain  rate  ea"  to  the  hold  stress  ah  for  stainless 
steels,  aluminum  alloys  and  ferritic  materials 


-  B'  sinh(/l<TA) 


(9.2) 


in  which  A,B'  and  n  are  material  parameters.  Miller  argues  that  during  a  creep  test  the  hardening 
variables  Y  and  R  assume  the  values  Yss  and  Ru ,  at  the  onset  of  steady  state  creep,  and  remain  con¬ 
stant  thereafter.  He  further  assumes  the  existence  of  a  function  such  that  during  steady  state 
creep 

°h-Yu^ 


4"  =  *t 

Substituting  expression  (9.3)  into  (9.2) 


eu"  *  B'  i  sinh 


=  <Jh 


Vh-Y* 


(9.3) 


(9.4) 


It  is  next  supposed  that  relation  (9.4),  which  holds  during  steady  state  creep,  is  also  appropriate  for 
non-steady  situations  as  well,  i.e., 

r 

,  .  i  a- Y 

All ) 


i"  «  AM  sinh 


(9.5) 


Further  recourse  to  empirical  argument  is  made  to  propose  that 


A * 


i 

Eq.  (9.1)  thus  acquires  the  specific  form 

t"  *  B'  Isinh 
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(9.6) 


a-Y 
R 

L  J 


3/2 


(9.7) 


for  the  inelastic  strain  rate. 


A  temperature  dependence  is  introduced  through  the  material  parameter  B  ,  based  upon  the 
observation  that  the  activation  energy  for  plastic  flow  varies  with  temperature.  It  is  expressed  by  the 
relations 

B'  =  Bo>,  (9.8) 


<j>  *  exp 


u  *  exp 


KQ 


d  >  0.6dm 


0<;O.6e„ 


(9.10) 


B  and  k  are  material  constants,  0m  is  the  melting  temperature  and  Q  is  the  activation  energy  neces¬ 
sary  for  the  plastic  deformation. 

The  hardening  variables  are  assumed  to  be  governed  by  the  following  rate  equations 

Y  -  H x'e"  -  *2(Y,G)  (9.11) 

R  «  Jf2|e"|  -*3(/?3,d) 

which  include  thermally  activated  softening  (or  hardening  recovery)  by  means  of  the  recovery  func¬ 
tions  +2  and  ♦3.  These  functions  are  taken  to  have  the  same  form  as  the  steady  state  creep  correla¬ 
tion  expression  (9.2).  Adding  the  elastic  strain  rate,  and  summarizing: 

e  -  «'  +e" 


3/21  1  « 


i"  -  Bo*d)|sinh  ^  J  j  |  Sg(T  -  Y) 

-  H{t”-HXBU8)  [sinh  [ax  |  Y  |  ]  ]  * Sg{Y) 

~  -  H2\i”\  -H2  flu(d)  [sinh(A2/?3)]" 
exp  — Q/kd\  6  2:  0.6dm 

«(«)  *  J  L  J 

exp  ([ -Q /0.6dm ]  [ln(-~p-—  +  1)])  d  JS  0.6dm 

0 


(9.12) 


[B,n,Hx,Hi,Ax,A2,Q,k,9m)  are  viscoplastic  and  thermal  material  parameters. 


10.  Valanis  Theory 
10a.  Intrinsic  Time  Scale 

The  derivation  of  the  constitutive  equations  of  Valanis’  ‘endochronic’  theory  of  viscoplasticity 
[53-61]  is  lengthy.  It  is  necessary,  however,  to  go  through  the  derivation  in  order  to  undersand  the 
theory’s  strengths  and  weaknesses.  Central  to  the  theory  is  the  notion  of  ‘intrinsic  time’  and  its  mea¬ 
sure,  which  is  developed  as  follows.  In  first  considering  rate  independent  plastic  deformation, 
Valanis  argues  “that  the  state  of  stress  in  the  neighborhood  of  a  point  in  a  plastic  material  depends 
on  the  set  of  all  previous  states  of  deformation  of  that  neighborhood,  but  it  does  not  depend  on  the 
rapidity  at  which  such  deformation  states  have  succeeded  one  another.  The  independence  of  stress  on 
the  rapidity  of  succession  of  deformation  states  is  achieved  by  introducing  a  time  scale  f,  called  the 
‘intrinsic  time’  which  is  independent  of  time  t,  the  external  time  measured  by  a  clock,  but  which  is 
intrinsically  related  to  the  deformation  of  the  material.”  This  relation  is  postulated  to  have  the  form 

df^dCPdC  (10.1) 


where  C  is  the  right  Green-Cauchy  deformation  tensor  and  P  is  a  fourth-order  material  tensor  which 
could  depend  on  C.  The  ‘intrinsic  time'  f  describes  the  length  of  the  path  traveled  by  the  material  in 
the  six-dimensional  deformation  space  and,  thus,  is  representative  of  the  deformation  history  of  the 
material. 

To  account  for  the  rate  dependence  exhibited  by  structural  materials  when  experiencing  non- 
recoverable  deformation,  a  second  time  scale,  L ,  the  ‘intrinsic  time  measure’,  coupling  the  intrinsic 
time  f  with  the  natural  time  t  is  introduced. 

dL2  •  df2 +g2dt2  (10.2) 

The  scalar  g  is  a  material  parameter.  Finally  an  ‘intrinsic  time  scale’,  z(L),  is  defined  such  that 

dz  *  'TTTTdL ,  ^>0,  OsL  <  oo  (10.3) 

f(L)  dL 

where  f(L)  is  another  material  function.  Since  P  and  /  are  material  functions,  the  intrinsic  time 
scale  z  is  viewed  as  representing  the  particular  inner  properties  of  the  material.  Hence  the  descriptive 
name  of  ‘endochronic’  given  to  the  theory. 

10b.  General  Theory 

The  general  form  of  the  theory,  which  allows  for  finite  deformation,  is  developed  from  contin¬ 
uum  thermodynamic  considerations  in  which  internal  variables  are  introduced  and  coupled  with  the 
Onsager-DeG root- Mansur  theory  of  non-equilibrium  thermodynamics.  Designating  t.ij  and 
♦  *  e— (hq  as  the  internal  energy,  entropy  and  free  energy  densities  per  unit  mass,  respectively,’"  q  as 
the  heat  conduction  vector  and  h  as  the  heat  source  density,  then  the  first  and  second  laws  of  thermo¬ 
dynamics  expressed  in  local  form  require  that  at  each  point  of  the  body  undergoing  a  thermo- 
mechanical  process 

pe  *  tr  (T  •  D)  -  div  q  +  p/r  (10.4) 


*In  this  section  the  symbol  t  is  used  to  designate  internal  energy  and  not  uniaxial  strain  as  in  other  parts  of  this  report. 


w  L-w¥i  nr;  vr  mr  jl"  itf  (  it  Vn  k' 


and 

tr(T- D)  -  pi'  -  p  if  9  -  -jq-  grad  0sO.  (10.5) 

d 

In  these  expressions  p  is  the  mass  density  at  a  point  of  the  current  configuration  of  the  body. 

For  elastic  material  response  it  is  assumed  that  the  material  behavior  can  be  specified  by  consti¬ 
tutive  equations  for  the  free  energy  ♦ ,  the  stress  T,  the  entropy  i\  and  the  heat  conduction  q,  when 
the  deformation  C  the  temperature  9  and  its  gradient  g  -  grad  9  are  presumed  to  be  known.  To 
describe  inelastic  response,  Valanis  argues  that  associated  with  non-recoverable  deformation  are 
numerous  processes  that  occur  at  the  atomic  scale  that  are  responsible  for  internal  entropy  generation. 
It  is  necessary,  therefore,  to  introduce  JV  additional  state  variables,  which  Valanis  refers  to  as  'inter¬ 
nal  hidden  variables’  (i.e.,  not  necessarily  observable)  in  order  to  completely  define  the  thermo¬ 
dynamic  state  of  the  material.  Thus,  according  to  this  view,  the  complete  thermodynamic  state  of  the 
material  requires  knowledge  of  C,  9,  g  and  V-many  unspecified  internal  variables,  /3,, /32, . . . ,  0,v ,  at 
each  point  of  the  material.  The  internal  variables  are  further  assumed  to  be  governed  by  N  differen¬ 
tial  equations  of  evolution. 

The  set  of  postulated  constitutive  equations  for  any  thermomechanical  process  is  thus  enlarged  to 
include  the  internal  variables  as  follows: 

♦  -  *(C,9,  . . ,0N) 

T  ■»  T(C,tf,g,)31 . fiN)  (10.6) 

V  *  ij(C,ff,g,0i,. . .  ,$y) 

q  *  q(C,9,g,0lt . 0N) 

and 

C,9,0u...,0n)  i  =  1,2,..., AT. 

Comment: 

1.  It  is  noted  that  Valanis  fails  to  include  the  temperature  gradient  g  as  a  variable  for  the  functions 
fi ,  even  though  it  is  included  in  all  of  the  other  constitutive  equations.  No  reasons  nor  justifica¬ 
tion  is  given  as  to  why  this  should  be  so. 

Following  the  analysis  of  Coleman  and  Guritn  [62]  substitution  of  Eqs.  (10.6)  into  the  first  and 
second  (aw  requirements  (10.4)  and  (10.5)  lead  to  the  following  conclusions:  The  constitutive  func¬ 
tions  for  ♦.Tandij  must  be  independent  of  the  temperature  gradient  g  as  an  explicit  variable.  In 


addition  the  free  energy  determines  the  stress  and  the  entropy  through  the  relations 

T  a,  2-£-— 

Po  3C  ’ 

d * 

where 

(10.7) 

♦  -  *(C,0,/3,,...  ,0*), 

and 

(10.8) 

a ♦  , 

d*  • 

(10.9) 

d(3,fi  “ 

— -pf20  i  *  1 . N(  no  sumoni). 

dpi 
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As  mentioned  above,  in  obtaining  expressions  (10.7)-(10.9)  Valanis  employs  the  arguments  of 
Coleman  and  Gurtin,  who  follow  the  Truesdell-Noll-Coleman  theory  of  non-equilibrium  thermo¬ 
dynamics.  However  in  developing  his  theory  further,  Valanis  next  turns  to  the  Onsager-DeGroot- 
Mansur  theory  when  he  makes  use  of  a  result  central  to  that  theory  of  non-equilibrium  thermodynam¬ 
ics.  According  to  the  inequality  (10.9)  the  product  of  (-  d  */d0,)  wrh  (d0,/dt)  results  in  an 
increase  of  the  entropy.  They  are  then  interpreted  as  ‘generalized  forces’  and  ‘generalized  fluxes’ 
when  Valanis  next  assumes  the  Onsager  relations  which  relate  the  generalized  forces  linearly  to  the 
generalized  fluxes,  that  is. 


3je_  ,di 3, 

30,  =  ij  dt  ’ 


i  =1 . N 


(  no  sum  on  /). 


(10.10) 


In  the  equations  the  b‘  can,  in  general,  be  functions  of  the  state  variables  C,0,  0, . 0N. 

When  viewed  as  a  set  of  differential  equations  for  the  unknown  internal  variables  0lt. . .  ,0N,  Eqs. 
(10.10)  are  non-linear  differential  equations.  A  linearization  is  introduced  by  further  assuming  the  b,j 
are  such  that 

b-j  =  k  b(J ,  (10.11) 


where  b,j  are  constants  and  k  >  0  is  a  parameter  that  relates  natural  time  t  to  the  intrinsic  time  scale 
z  whereby 


dz 


(10.12) 


Insertion  of  Eqs.  (10.11)  and  (10.12)  into  Eq.  (10.10)  then  produces  the  following  set  of  linear  dif¬ 
ferential  equations  for  the  internal  variables  0j : 


3+ 

30, 


+  h.-, 


d0j 

dz 


■  0, 


1,2,. 


N  {  no  sum  on  j). 


(10.13) 


The  system  of  linear  differential  Eqs.  (10.13)  for  determination  of  the  0j  is  set  in  terms  of  the  intrin¬ 
sic  time  scale  z,  and  replaces  the  system  of  generally  non-linear  differential  equations  (10.6);  set  in 
terms  of  the  natural  time  t . 

% 

With  equations  (10.13),  the  constitutive  description  of  the  entropy  producing  (dissipative) 
material  behavior  can  now,  in  principle,  be  considered  to  be  complete.  One  needs  only  to  expand  the 
free  energy  ♦  as  a  Taylor  series,  which,  in  principle,  allows  solution  of  the  system  of  equations 
(10.13)  for  the  internal  variables  0,,... , 0 These  values  when  substituted  into  Eqs.  (10.6)2  g've 
the  desired  stress-deformation  relations. 


In  summary: 
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3* 

30, 


0, 


(10.14) 


j  =  1,2,. ..  ,1V  (no  sum  on  j). 


Comments: 


2. 


If  in  the  constitutive  equations  (10.6)s  the  temperature  gradient  g  =  grad  6  had  been  included  as 
an  independent  variable  as  it  was  for  the  other  constitutive  equations,  then  in  place  of  the  dissi¬ 
pative  inequality  (10.9)  there  would  be 


-  —  +  — qq  >  0 

d0j  ‘  p# 


instead.  It  appears  in  light  of  the  Onsager  type  relations  (10.10),  which  are  suggested  by  the 
simpler  inequality  (10.9),  and  which  make  possible  the  expressions  (10.13),  that  herein  lies  the 
rationale  for  Valanis’  omission  of  g  =  grad  6  as  a  possible  variable  for  the  rate  functions  /,  at 
the  outset.  The  above  inequality  would  preclude  introduction  of  the  Onsager  like  relations 
(10.10). 


3.  A  unique  solution  to  the  set  of  N  first  order  differential  equations  (10.6)s  or  (10.13)  for  the  N 
internal  variables  j3t, . . .  ,0N,  requires  stipulation  of  a  set  of  N  initial  conditions  for  these  vari¬ 
ables.  This  presents  some  difficulty  in  light  of  the  interpretation  of  these  variables  as  “hidden” 
unspecified  variables.  How  can  physically  relevant  initial  conditions  be  assigned  on  to  set  of 
variables  if  they  are  not  defined? 


10c.  Explicit  Constitutive  Equations  for  Small  Deformations 


The  hidden  variables  heretofore  undefined  as  to  generic  nature,  are  now  assumed  to  be  N 
second-order  symmetric  tensors  0-J  =  0J},  i  ,j  =  1, 2, 3,  a  =  1, . . . ,  N  (we  hereafter  refer  to  a  rec¬ 
tangular  coordinate  system  and  resort  to  a  tensor  component  notation).  For  small  deformation  and 
uniform  temperature,  Eqs.  (10.14),  after  appropriate  modification4',  simplify  to 

*«*(  Ey,0'y,...,0#)  (10.15) 


a* 

d/3 


dZ, 


biju 


£0S_ 

dz 


-  0 


a  =  1 , . . . ,  N( a  not  summed  ). 


(10.16) 


(10.17) 


The  components  b^  are  the  components  of  the  N-many  fourth-order  tensors  b“,  a  -  1 . N , 

which  are  considered  to  be  material  tensors. 


The  free  energy  is  expanded  as  a  Taylor  series  up  to  terms  of  the  second  order.  (The  terms  of 
the  first-order  vanish  to  satisfy  the  zero  deformation  of  the  initial  undeformed  state). 


2  ^ijkl  EjjZ, 


kl 


L  BfaZijPS 

am\ 


J  Z  C-ju  0fj  05 . 

Om  1 


(10.18) 


Here  A,  Ba,  C°  are  material  tensors.  Assuming  material  isotropy,  initially  and  throughout  any  subse¬ 
quent  deformation,  the  fourth-order  material  tensors  are,  therefore,  isotropic  tensors  and  can  therefore 
be  given  the  representation 

AjjU  *  A 1 5y  bu  +  A  2  &n  bji 

Bfa  •B?6,J6u+B?&ik8JI  (10.19) 

CijU  *  C{*  bijbki  +  Cf  bikbjt 
btjki  =  Sjj&id  +  bf  bik6jt 

'Having  the  internal  variables  as  second-order  tensor  quantities,  it  is  correspondingly  necessary  the  the  N  -nany  constant 
coefficients  in  Eqs.  (10.13)  be  replaced  by  components  B fa  of  the  N-many  fourth-order  tensors  B“. 


where  A\,Ai,Bf  ,Bg  ,Cf  ,Cg ,bf  ,bf ,  a  =  t N  are  material  constants  and  5jk  are  the 
Kronecker  delta  symbols.  Thus 


♦  =-  +  JE  (fl?5iy5u  +  Bfbik 6#)% ft 

+  J  J,  (C"M«  +  C2“6,y57V)  0g0g  , 


(10.20) 


which  gives 


=  5f  +  B?Zij  +  Cftyi 3&  +  Cf0fj. 

Wij 


(no  sum  on  a).  (10.21) 


Also,  from  (10.18)  and  (10.19)4 


d0S  d0&  d0g 

a  **  _  i.  or  w  c  i  f_  a 

~  bi~8ij  +b*^r 


With  these  results,  Eq.  (10.14)4  now  takes  the  form 


(  no  sum  on  a) . 


(10.22) 


d03t  d0i j} 

+  bf- f-  5  y  +  B|-Ey  +  Cftf  +  b?-±  =  0, 


(10.23) 


for  each  value  of  a  =  1,2,...  ,1V.  Solution  of  this  set  of  first-order  differential  equations  is  facili¬ 
tated  by  a  decomposition  into  spherical  and  deviatoric  parts  as  follows: 


d0& 

BgL*  +  CM  +  &0“-~  =  0, 

O 

°  °  d0g 

BfLij+CfB*  +  bf-£-  =  0, 


where 


K  -  Bf]  ~  yfl&fy. 


Bg  =  3fl?  +  Bf,  eg  =3 Cf  +  C2“,  60“  =  36“  +  62“, 
The  general  solutions  of  Eqs.  (10.24)  and  (10.25)  are  given  by 


D0  rz 


"2  <■ 

W  -  -  tK 


Ettfe'Vfc'  +  Dae 


(10.24) 


(10.25) 


(10.26) 


(10.27) 


(10.28) 


X,j(z’)dz'  +  Fae 


(10.29) 


r«,i 


in  which  Da  and  Ea,  a  =  1, . . .  ,  N,  are  arbitrary  constants  and 

K  =  Cg/bg,  Pa=Cg/bg. 


(10.30) 


Comment: 

4.  Since  no  initial  conditions  for  the  set  of  differential  equations  (10.24)  and  (10.25)  are  given 

(most  probably  because  the  internal  tensor  variables  . 0N  are  not  defined,  thus  there  is 

no  physical  motivation  available  to  allow  assumption  of  an  appropriate  set)  the  set  of  arbitrary 
constants  Da  and  Ea  cannot  be  determined.  In  as  much  as  the  equations  that  follow  do  not 
include  the  Dae  and  Fae  P“~  parts  of  the  solutions  (10.28)  and  (10.29)  one  can  only 
conclude  that  the  arbitrary  constants  Da  and  Fa  were  all  set  to  zero,  although  no  argument  or 
explanation  as  to  why  this  should  be  so  is  given. 


Eqs.  (10.28)  and  (10.29)  (with  the  terms  including  Da  and  Fa  hereafter  omitted)  determine  the 
undefined  internal  variables  in  terms  of  the  strain  and  the  intrinsic  time  scale.  Upon  substitution  into 
the  free  energy  form  (10.18),  calculation  of  stress-strain  relations  (10.16)  may  then  follow.  From 
(10.18)  and  (10.16) 


T(j  —  A^L^bij  +  A2'E‘ij  +  £  \bi0Sc  &ij  +  Bg  0 “  . 

at*  1  L  J 


which  can  be  decomposed  into  the  deviatoric  and  spherical  parts 


Tij  =  +  Bg  kj 


(10.31) 


(10.32) 


Tkk  *  AgE/dc  +Bg  Bkk, 

A0  m  3Al+A2,  Bg  =  3fl“  +  B2 


(10.33) 

(10.34) 


Substitution  of  (10.29)  into  (10.32),  followed  by  integration  yields 


»  BgBg  o  ,  BgBg  dZH 

T(j  =  (/4  +  f  -=-=-e  P“U  z)— -ftfc'. 

v  2  Cg  '  ,J  Jo  eg  dz' 


(10.35) 


Designating  by  H(z)  the  unit  step  function  and  defining 


BgBg  n  BgBg  . 

Mz)  =  a2— 2^2-  H(z,  +  E  -2^2-  «  " 

Cg  0-1  C g 


(10.36) 


Eq.  (10.35)  takes  the  form 


•  dL,°. 

T‘i  • 


(10.37) 


in  which  n(z)  is  interpreted  as  a  shear  modulus.  By  a  similar  calculation  in  which  Eq.  (10.28)  is  sub¬ 
stituted  into  Eq.  (10.33),  integrated,  and  where  the  definition 


"  BgBg  n  BgBg  _x. 

K(z)  =  M0  -  E  — —  H(z)  +  E  -2-2-  e 


eg 


■'  eg 


(10.38) 


(10.39) 
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is  introduced,  the  spherical  part  of  the  stress  tensor  takes  the  form 

y  n»  -£KU-z')-£rZudz’, 

in  which  K(z)  is  considered  to  be  a  bulk  modulus. 

After  a  re-labeling  of  the  expressions  for  the  shear  and  bulk  moduli,  in  summary: 


O 

dZ, 


dZ 


T^^z-z^dz'^^Kiz-z^dz' 


dz 


dz' 


N 


n(z)  =  H 0+  Z  nae 

a”  1 

N 


K(z)  =  K0+  Z  Kae 

at*  l 


(10.40) 


{/lo,  K0,  fia,  Ka,  pa,  \a,  a  =  1 , . . .  JV  )  are  material  parameters 
z  »  intrinsic  time  scale  determined  by  Eqs.  (10. 1 )-( 10.3). 


The  theory  is  carried  further  for  rate- independent  deformation  by  assuming  the  parameter  g  in 
Eq.  (10.2)  is  zero.  The  intrinsic  time  measure  is  then  (for  small  deformation) 


dz2  =  dZPdZ  =  PijudZijdZu. 


(10.41) 


The  material  tensor  P  is  isotropic  [cf.  Eqs.  (10.19)] 

Pijki  *  kAj&kl  +  Mi*5,/ 


so  that 


L  =Jo"  k,(dEtt)2  +  kJLijdL'j 


1/2 


(10.42) 


The  material  function  /(z)  of  Eq.  (10.3)  is  taken  as 

f(z)  =  1  +  0  L  0>O 


(10.43) 


from  which  it  follows  that 


z  =  ~-ln(l  +0L). 


(10.44) 


The  constants  k2  and  0  are  additional  material  parameters. 


The  strain  input  into  Eq.  (10.42)  determines  the  intrinsic  time  measure  L,  which  then,  through 
Eq.  (10.44),  specifies  the  intrinsic  time  scale  z. 
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Comments: 


For  inelastic  and  plastic  deformations  that  exhibit  anisotropic  hardening  e.g.,  kinematic  or  direc¬ 
tional  hardening,  the  theory  given  by  Eqs.  (10.40)  is  clearly  inappropriate. 

The  undefined  internal  variables  a  =  1,2,. . . ,  N  are  solved  for  and  expressed  in  terms  of 
six  sets  of  N-many  material  parameters  a  =  I,...,  N  and  the 

intrinsic  time  scale  z.  These  are  eventually  lumped  together  as  the  material  functions  (10.40)2  3 
involving  4N+2  undefined  parameters,  where  the  number  N  is  left  open.  Since  the  internal 
variables  themselves  are  not  identified,  the  question  of  what  should  be  the  value  of  N  is  only 
vaguely  alluded  to.  In  recent  applications  of  the  theory  the  number  has  ranged  from  /V  =  1  to 
N  =  oo.  What  does  it  mean  to  say  that  in  one  situation  one  internal  variable  is  necessary  to 
describe  the  material  behavior  where  as  in  another  an  infinite  number  of  them  are  necessary? 
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COMPARISONS  BETWEEN  THE  MODIFIED  CHABOCHE  AND  THE 
BODNER  THEORIES  FOR  A  HIGH  TEMPERATURE  APPLICATION 


h 


i 


* 


> 


11.  Introduction 

The  modified  Chaboche  and  the  Bodner-Partom-Stouffer  theories  were  selected  to  model  the 
behavior  of  the  alloy  Inconel  at  a  temperature  of  1200°F.  Each  theory  appears  to  be  the  best 
representative  from  the  category  of  theories  in  which  they  were  grouped.  Specifically,  the  behavioral 
characteristics  of  Inconel  at  1200°F,  as  predicted  by  both  theories,  are  examined  for  uniaxial  mono¬ 
tonic  tensile  loading  at  varying  strain  rates,  for  sudden  increase  and  decrease  of  strain  rate  in  the  ine¬ 
lastic  range,  for  creep  and  stress  relaxation  behavior,  for  load-unload-reload  at  varying  strain  rates, 
and  for  cyclic  loading  [34]. 


12.  Viscoplastic  Material  Parameters  for  Inconel  at  1200 °F 


For  uniaxial  loading  the  modified  Chaboche  equations,  which  consist  of  Eqs  (5.14)  and  (5.19), 
are  repeated  here  in  combined  form. 

<  ■  «'  +«" 


K 


Sg(a-Y)  F  >0 
F  =  <;  0 


F  -  | a—Y\  -R(pe)  Sg  (a-Y)  = 


a-Y 

\a-Y\ 


=  ±1 


(12.1) 


R(p,'e)  =  q+  /?0(e) -q  1  p  \dt"\ 

L  J  o 

^  =c(ak"-Y \i”\  )  -  a  |  Y\m  Sg(Y) 

[K,n,b,q,c,a,y,m\  are  viscoplastic  material  parameters. 
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For  high  temperature  application  in  which  hardening  recovery  (thermal  softening)  is  expected,  all 
eight  of  the  viscopiastic  material  parameters  require  determination.  In  addition,  the  form  of  the 
dependence  of  the  initial  yield  stress  on  the  strain  rate,  /?0(«),  must  also  be  ascertained. 

At  1200°F,  the  available  experimental  data  for  Inconel  [35-38]  shows  that  the  alloy  exhibits 
strain  rate  sensitivity  only  over  the  small  strain  rate  range  «  <  10-4  sec  The  material  also  shows 
typical  creep  response  (primary,  secondary,  tertiary  creep  stages),  and  cyclic  softening  for  approxi¬ 
mately  10  to  20  percent  of  fatigue  life  before  cyclic  stabilization.  The  stress-strain  response  tends  to 
a  constant  stress  value  at  1.0  to  1.5  %  of  strain. 

The  manner  by  which  values  for  the  parameters  were  evaluated  from  the  data  of  monotonic 
loading  tests,  cyclic  loading  and  creep  tests,  is  described  in  References  [33,63],  with  the  results  tabu¬ 
lated  below. 


Parameter 

Description 

Value 

E 

Young’s  modulus 

24.73  x  lfr*  ksi 

V 

Poisson’s  ratio 

0.336 

a 

Saturation  value  of  kinematic 
hardening  variable 

30.00  ksi 

c 

Kinematic  hardening  exponent 

350.0 

7 

Coefficient  of  hardening  recovery 

0.4  x  10-,° 

m 

Hardening  recovery  exponent 

7.00 

K 

Overstress  parameter 

155.0  ksi  sec1''" 

n 

Strain  rate  sensitivity  parameter 

5.10 

q 

Saturation  value  of  isotropic 
hardening  variable 

50.00  ksi 

b 

Isotropic  hardening  exponent 

3.75 
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At  the  elevated  1200°F  temperature,  the  experimentally  observed  strain  rate  sensitivity  of  the 
initial  yield  stress  was  primarily  over  the  low  strain  rate  range  10~9  <  e  <  10~4  sec  ~ 1 ,  where  the 
data  is  well  correlated  by  Eq.  (15.21),  [32,63] 


*o(<)  =  ^00^  1  +(1+  exp  £  a,  One)' 

I  ,-o 


where 


lim  R(p ,()  —  R0(Q)  —  /?oo 

(p.t)— o 


(12.2) 


represents  the  yield  stress  at  quasi-static  strain  rate  and  the  a,  are  non-linear  regression  coefficients. 

The  procedures  for  determination  of  the  material  parameters  appearing  in  the  Bodner-Partom- 
Stouffer  viscoplastic  constitutive  equations: 

e  =  e'  +«" 


*’  “  j* 


«"  =  ~77Do  exp  -  —r 


z  -  zo  +  q  I0  Z(r)dr  +  (1  -  q)Sg(k")  fetf-C-dr 


(12.3) 


Z  =  m(z, -Z)W"  ~  Az\ 


W  =  [  af"dt  =  J*  adt" 
J0  '0 


DQ,n,Zo,Z\,Zi,m,A,r,q  >•  are  viscoplastic  material  parameters 


are  detailed  in  Ref.  [38]  and  the  values  are  shown  below. 


Parameter 

Description 

Value 

E 

Young's  modulus 

23.55  x  10J  ksi 

V 

Poisson’s  ratio 

0.336 

m 

Rate  of  work  hardening 

2.9 

2o 

Initial  value  of  hardness 

235.1  ksi 

Z| 

Saturation  value  of  hardness 

260.2  ksi 

Z2 

Minimum  recovery  value  of 
hardness 

104.1  ksi 

A 

Hardening  recovery  parameter 

1.5  x  10-3 

r 

Hardening  recovery  exponent 

7.0 

n 

Strain  rate  sensitivity  parameter 

3.0 

Do 

Limit  value  of  the  inelastic  strain 
rate  (assigned  value) 

106  sec-1 

q 

Hardening  ratio  (assigned  value) 

-0.05  to  -0. 10 

13.  Experimental  and  Qualitative  Comparisons 
13a.  Experimental  Comparisons 


Several  comparisons  of  the  predictive  capability  of  the  two  theories  with  experimental  data  (to 
the  extent  permitted  by  the  availability  of  the  data)  are  shown  in  Figures  13.1  thru  13.3  for  monotonic 
tension,  for  a  fully  reversed  load  cycle  and  for  minimum  or  secondary  creep  rate.  The  agreements 
shown  between  both  sets  of  predictions  and  the  experimental  results  are  reasonably  good. 

It  is  of  interest  to  observe  that  the  Bodner  theory  which  does  not  employ  a  yield  condition,  and 
which  allows  for  inelastic  deformation  at  all  levels  of  stress,  nevertheless  shows  a  stress-strain  curve 
with  a  rather  sharply  defined  transition  from  linear  (and  presumably  reversible  elastic)  behavior  to 
non-linear  behavior,  at  a  stress  level  close  to  what  would  be  considered  as  a  yield  stress  in  those 
theories  that  assume  the  existence  of  a  yield  condition. 

13b.  Strain  Rate  Effects 

The  modification  introduced  into  the  Chaboche  theory  allows  the  theory  to  be  receptive  to  the 
effect  of  strain  rate  on  the  initial  yield  stress.  Figures  13.4  and  13.5  illustrate  the  effect  of  strain  rate 
on  the  stress-strain  response  of  Inconel  at  the  elevated  temperature  as  predicted  by  the  two  theories 
for  three  different  constant  strain  rates:  5x  10~5,  5x  10-7,  and  5x  10-9  sec'1.  Both  theories 
predict  a  pronounced  strain  rate  effect  in  the  strain  range  considered.  In  both  cases  the  stress 
saturates  at  relatively  low  strains,  which  is  consistent  with  what  is  observed  experimentally  for 
superalloys  at  high  temperature.  Also,  the  inelastic  parts  of  the  stress-strain  curves  tend  to  approach 
one  another  as  the  strain  rates  increase  in  constant  ratio,  a  result  that  is  consistent  with  experimental 
observation  in  general. 

A  qualitatively  similar  stress-rate  effect  can  also  be  simulated  by  both  theories.  This  is  shown 
by  the  stress-strain  responses  obtained  under  stress  controlled  conditions.  The  stress  rates  used  were 
E  (elastic  modulus)  times  the  corresponding  strain  rates.  The  results  show  the  stress-controlled 
predicted  curves  to  lie  entirely  above  the  corresponding  strain-controlled  curves  throughout  the  non¬ 
linear  range.  This  difference  in  the  predicted  responses  is  attributed  to  the  non-linearity  of  the  stress- 
strain  rate  constitutive  relations. 

The  relative  strength  of  the  strain  rate  effect,  that  is  the  change  in  the  stress-strain  response  as 
the  strain  rate  is  increased,  appears  to  be  greater  for  the  Chaboche  theory  compared  to  the  Bodner 
theory  as  evidenced  by  the  relative  stress  difference  between  the  stress-strain  curves  for  the  e3  and  e, 
strain  rates. 

13c.  Strain  Rate  History  Effect 

The  response  of  viscoplastic  materials  can  be  influenced  by  the  strain  rate  history  of  the 
material.  This  can  be  shown  by  strain  rate  jump  tests.  Experiments  of  this  kind  on  polycrystalline 
materials  generally  show  the  following  characteristics:  As  the  strain  rate  from  a  plastic  state  is  sud¬ 
denly  increased  the  initial  response  of  the  material  is  predominately  elastic,  followed  by  yield  and 
subsequent  hardening  during  which  the  jump  response  approached  asymptotically  the  monotonic 
response  at  the  higher  strain  rate.  Whether  the  jump  response  approaches  the  monotonic  response 
from  the  low  side  or  the  high  side,  has  been  shown  to  depend  upon  the  ratio  of  the  strain  rate  before 
and  after  the  jump,  and  on  the  crystalline  structure  of  the  material.  For  example  for  FCC  and  HCP 
materials  the  jump  response  is  lower  than  the  monotonic  response,  while  for  BCC  materials  the 
reverse  takes  place. 


STRESS,  ksi 


STRAIN,  % 


R«.  13.1  —  Experimental  and  predicted  tensile  stress-strain  curves  at «  *  1.333  x  10~5  sec' 
using  the  modified  Chaboche  and  Bodner  theories. 


-1.2  -0.8  -0.4  0.0  0.4  0.8 


STRAIN,  % 

Fig.  13.2  —  Experimental  and  predicted  Strain-controlled  fully-reversed  cyclic 
behavior  at  <  •  4  x  10"’  sec ' '  and  -  2%  using  the  modified  Chaboche  theory; 
expert  mental  peak  stress  only  is  available. 
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/  cjsSxlO"*5  sec-1 

“/ STRAIN  CONTROL  e2  =  5  x  1 0  ~  7  sec  " 1 
/  63  =  5x10~9  sec “ 1 


/  I  or-j  —  1 .237  x  10°  ksi/sec 

f— STRESS  CONTROL  la2  —  1  -237  x  1 0~2  ksi/sec 

I  a3  =  1 .237  x  1 0~4  ksi/sec 


0.00  0.30  0.60  0.90  1.20  1.50  1.80 

STRAIN,  % 

Fig.  13.4  —  Predicted  tensile  response  at  different  rates  employing  the 
modified  Chaboche  theory. 
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Fig.  13.5  —  Predicted  tensile  response  at  different  rates  employing 
Bodner's  theory. 


Figures  13.6  and  13.7  illustrate  the  strain  rate  jump  response  predicted  by  the  two  theories  as 
the  strain  rates  are  increased  (decreased)  by  one  order  of  magnitude,  from  10-6  sec-1  to  10~5  sec-1. 
Both  theories  indicate  an  initial  elastic  response  immediately  following  the  jump  in  the  strain  rate,  fol¬ 
lowed  by  hardening  and  a  gradual  approach  to  the  monotonic  stress-strain  curve.  Both  demonstrate 
the  strain  rate  history  effect  in  the  sense  that  the  jump  response  does  not  immediately  approach  the 
monotonic  response  because  of  its  prior  inelastic  pre-strain  history  at  the  lower  (upper)  strain  rate. 
The  approach  of  the  jump  response  to  the  monotonic  response  occurs  more  rapidly  according  to  the 
Bodner  theory . 

13d.  Creep  Behavior 

There  is  experimental  data  that  indicates  that  the  creep  behavior  of  structural  materials  can  be 
significantly  influenced  by  the  rate  of  stress  that  is  employed  in  reaching  the  hold  stress  ah  in  a  uniax¬ 
ial  creep  test  [64].  Figures  13.8  and  13.9  demonstrate  that  both  viscoplastic  theories  can  model  this 
observed  effect,  where  the  higher  stress  rate  (ax  =  20  ksi/sec)  to  the  hold  stress  shows  considerably 
greater  primary  creep  strain  than  the  lower  stress  rate  (a2  =  20  ksi/sec),  prior  to  the  attainment  of  the 
minimum  creep  rate  (secondary  creep).  The  minimum  creep  rate  appears  to  be  unaffected  by  the 
stress  rate  to  the  hold  stress. 

This  effect  may  be  explained  in  the  following  manner.  The  yield  stress  for  the  smaller  stress 
rate  is  lower  than  for  the  higher  stress  rate,  i.e.,  Rq  (e2)  <  Rq  («()  for  a2  <  ai-  The  hold  stress  at 
ah  =  140  ksi  is  well  above  both  yield  stress  values.  Yield  therefore  initiates  well  before  the  hold 
stress  level  is  attained,  and  by  the  time  the  creep  hold  stress  is  reached  the  material  has  hardened. 
For  the  smaller  stress  rate  with  the  correspondingly  lower  yield  stress,  the  material  will  experience  a 
greater  total  strain  in  getting  to  the  hold  stress  than  for  the  higher  stress  rate.  However  it  will  also 
have  experienced  greater  hardening.  Therefore  when  creep  begins  the  material  with  the  greater  prior 
hardening  will  undergo  creep  at  a  smaller  rate  and  thus,  for  a  given  time  interval,  have  a  smaller 
creep  strain,  even  though  the  total  strain  is  greater. 

Figure  13.10  shows  the  predicted  creep  response  after  5  1/2  hours  at  a  hold  stress  of  125  ksi. 
The  Bodner  creep  curve  shows  primary  and  secondary  creep  stages  although  the  primary  creep 
predicted  is  too  small  to  appear  on  the  scale  of  this  plot.  The  Chaboche  creep  curve  shows  the  con¬ 
ventional  primary  and  secondary  stages  of  creep.  However  it  also  shows  gradually  increasing  creep 
rates  prior  to  the  tertiary  creep  stage. 

Numerical  analysis  [63]  indicates  that  with  the  initiation  of  creep,  at  1200°F  Inconel  starts  to 
harden  kinematically  and  to  soften  isotropically,  where  the  rate  of  hardening  is  considerably  greater 
than  the  rate  of  softening,  Y  >  |  —R  \ .  This  is  so  because  R  depends  upon  the  inelastic  strain  accu¬ 
mulation  which,  during  the  initial  stages  of  the  creep,  is  very  small.  Thus  with  overall  hardening 
there  will  be  decreasing  creep  rate,  which  is  the  primary  stage  of  creep.  During  the  secondary  creep 
stage,  where  the  creep  strain  rate  is  constant,  it  can  be  easily  seen  from  Eq.  (5.17)  that  Y  =  |  -R  |  . 
The  rate  of  kinematic  hardening  is  equal  to  the  rate  of  isotropic  softening  while  the  creep  strain  rate  is 
constant.  Thus,  when  the  creep  strain  rate  starts  increasing  gradually  before  the  beginning  of  tertiary 
creep,  this  must  apparently  be  caused  by  material  softening,  that  is,  Y  <  \  -R  | ,  the  rate  of  isotropic 
softening  now  exceeds  the  rate  of  kinematic  hardening. 

Instances  of  gradually  increasing  creep  rates  following  a  period  of  steady  state  creep,  prior  to 
the  tertiary  stage  of  rapidly  increasing  creep  rates  associated  with  material  damage,  have  been 
observed  experimentally  for  structural  materials  at  elevated  temperatures  [65,66], 
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3.9  —  Predicted  creep  responses  at  140  ksi  and  different  stress  rates 
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Fig.  13.10  —  Predicted  creep  response  at  125  ksi  by  the  modified 


13e.  Relaxation  Behavior 


Structural  alloys  at  high  temperature  show  some  degree  of  stress  relaxation  at  fixed  strain,  with 
the  stress  normally  decreasing  at  a  decreasing  rate  until  a  stabilized  value  is  reached.  At  a  constant 
value  of  the  total  strain,  e  =  e'  +e",  it  follows  that 


e  = 


(13.1) 


so  that 


a  =  -  E'i” 


(13.2) 


which  indicates  that  the  stress  is  decreasing,  that  is,  relaxing.  At  stabilization  b  =  0. 

For  the  Chaboche  theory,  by  means  of  Eqs.  (12.1)3  and  (13.2)  the  condition  for  stabilization 
requires  that 


f  (a-Y)-R 

l  * 


=  0, 


(13.3) 


from  which 


a  =  Y  +  R 


(13.4) 


as  the  value  for  the  stress  when  the  stress  relaxation  has  ceased.  For  Bodner’s  theory  the  condition 
b  —  0  leads  to 


D0  exp 


=  0, 


(13.5) 


which  further  requires  that 


* 

n  + 1 

r  *v 

z 

2/i 

a 

(13.6) 


This  condition  cannot  be  satisfied  at  fixed  strain  unless  u— 0.  In  other  words,  the  stress  continues  to 
relax  with  time  until  it  vanishes  according  to  Bodner's  theory. 

Figure  13.11  shows  the  stress  relaxation  prediction  for  each  theory  for  a  loading  history  which 
strains  the  material  at  5  x  10"5  sec"1  until  1.0%  strain  is  attained  and  thereafter  held  fixed  for  a 
period  of  4  hours.  Both  theories  predict  a  plausible  early  response.  However,  the  Chaboche  stress 
relaxation  curve  shows  a  maximum  stress  reduction  of  approximately  15%,  with  stabilization  occur¬ 
ring  at  about  2  hours  after  the  strain  was  held  constant,  whereas  the  Bodner  stress  relaxation  curve,  as 
indicated  above,  continually  decreases  to  zero. 


13f.  Load-Unload  Reload  Behavior 

For  rate-independent  plasticity,  unloading  from  a  plastic  state  is  invariably  elastic.  This  is  not 
necessarily  the  situation  with  viscoplasticity.  The  initial  part  of  unloading  may  be  inelastic,  depending 
upon  the  ratio  of  the  loading  strain  (or  stress)  rates  to  the  unloading  strain  (or  stress)  rates.  Figures 
13.12  and  13.13  illustrate  the  predicted  response  of  the  two  theories  at  different  unloading  strain  rates 
after  a  strain  of  0.75%  was  reached.  The  symbol  5  signifies  the  ratio  of  the  loading  rate  to  the 
unloading  rate.  Thus  for  a  loading  strain  rate  of  5  x  10“7  sec-1,  the  indicated  5  values  correspond 
to  the  following  unloading  strain  rates: 

5,  =  -102-  e,  =  -  5  X  10“5jec-1 

62  =  —  1  —  C2=-5x  \Q~n  sec  “ 1 

63  =  -10"2  —  €3  =  -5x  10"9  sec~l. 

Qualitatively  speaking,  both  theories  predict  the  same  behavior.  They  show  that  the  initial 
unloading  behavior  is  rate  dependent.  For  63  =  — 10~2,  the  unload  strain  rate  is  two  orders  of  magni¬ 
tude  smaller  than  the  load  strain  rate,  and  the  initial  portion  of  the  unload  deformation  is  inelastic. 
Along  the  vertical  portion  of  the  unload  curve  the  total  strain  is  constant,  so  that  the  inelastic  strain  is 
increasing  in  magnitude  equal  to  the  decrease  in  the  elastic  strain,  e"  =  —  e'.  This  also  means  that 
the  material  is  hardening  during  this  part  of  its  unloading  history,  before  passing  into  elastic  unload¬ 
ing. 

Figures  13.14  and  13.15  give  the  predicted  responses  for  loading  and  unloading  at  different  rates 
of  stress.  The  stress  rates  used  were  determined  by  multiplying  the  strain  rates  that  produced  Figs. 
13.12  and  13.13  by  the  elastic  modulus.  The  loading  stress  rate  is  a  -  1.24  x  10~2  ksi/sec.  Thus 
the  indicated  a  values  correspond  to  the  following  unloading  stress  rates: 

71  =  102  —  er,  =  —  1.24  ksi/sec  . 

72  =  -1  —  ff2  =  -  1-24  x  10“2  ksi/sec. 

73  =  -10“2  —  =  -  1.24  x  10“4  ksi/sec. 

Both  theories  show  that  under  stress  control  the  initial  unloading  behavior  is,  as  to  be  expected,  also 
rate  dependent.  At  73  =  -  10“ 2  the  initial  unload  response  is  very  nearly  akin  to  creep  behavior.  A 
large  inelastic  strain  increase  accompanies  a  small  reduction  of  stress  before  the  elastic  unloading 
begins  to  take  place. 


0.00  0.15  0.30  0.45  0.60  0.75  0.90 

STRAIN,  % 


F'g-  13.12  —  Predicted  loading-unloading  response  at  different  unloading  strain  rates 
by  the  modified  Chaboche  theory. 
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Figures  13.16  and  13.17  indicate  qualitatively  similar  load-unload- reload  viscoplastic  behavior, 
and  illustrate  the  effect  of  the  unloading  rate  on  the  subsequent  reloading.  Shown  are:  (a)  A  load- 
unload-reload  cycle  at  the  same  absolute  strain  rate  | e (  =  5  x  10-2  sec  and  (b)  loading  and 
reloading  at  the  same  strain  rate  with  unloading  at  a  smaller  strain  rate.  In  both  cases,  the  strain  rate 
for  loading  was  chosen  at  5  x  10-5  sec-1  with  an  unloading  rate  of  -5  x  10~7  sec-1  selected  for 
the  second  case.  In  the  first  case,  5  =  -1  at  unloading,  the  initial  unloading  response  is  inelastic, 
approaching  very  rapidly,  however,  elastic  unloading.  The  reloading  response  correspondingly  does 
not  exhibit  any  discernible  difference  when  compared  to  the  monotonic  reloading  at  the  same  loading 
rate.  In  the  second,  6  =  — 10-2  at  unloading,  significant  inelastic  deformation  occurs  initially  upon 
unloading  over  the  vertical  portion  of  the  unload  curve.  Consequently,  because  of  the  accompanying 
strain  hardening,  the  subsequent  reloading  is  considerably  different  from  the  monotonic  response  at 
the  original  loading  rate. 

13g.  Cyclic  Load  Behavior 

At  1200  °F  the  cyclic  behavior  of  Inconel  shows  that  is  softens  for  approximately  10%  to  20% 
of  the  fatigue  life  prior  to  a  period  of  cyclic  stabilization.  In  the  Chaboche  theory  cyclic  softening 
and  subsequent  stabilization  are  modeled  by  means  of  the  isotropic  hardening  variable  R,  with  its 
saturation  q  below  its  initial  yield  value  R0(e).  For  Bodner’s  theory  the  parameter  q,  which  propor¬ 
tions  the  hardening  into  isotropic  and  directional  parts  must,  because  of  the  softening,  be  assigned  a 
negative  value.  It  was  estimated  in  Ref  [38]  that  for  Inconel  at  1200°F  a  value  of  q  between  -0.05 
and  -0. 10  is  appropriate. 


The  first  few  hystersis  loops  of  a  strain-controlled,  fully  reversed  cyclic  test  at 
|  e  |  =  4  x  10-5sec'1  strain  rate  and  At  =  2.0%  strain  range  as  predicted  by  the  two  theories  are 
shown  in  Figs.  13.18  and  13.19.  These  curves  show  that  the  Chaboche  theory  simulates  the  actual 
softening  quite  realistically,  as  indicated  by  the  continuous  decrease  in  the  stress  range  and  a 
corresponding  increase  in  the  inelastic  strain  range.  The  predictions  shown  by  the  Bodner  theory,  on 
the  other  hand,  are  not  typical  of  cyclic  softening.  The  predicted  hystersis  loops  have  almost  constant 


stress  and  strain  amplitudes.  The  predicted  curves  of  Fig.  13. 19  are  for  a  value  of  q 
values,  ranging  form  -0.05  to  -0.20,  were  also  tried  with  less  favorable  results. 


=  -  0. 10.  Other 


Parameter  Description _ 

E  Young's  modulus 

v  Poisson’s  ratio 


Saturation  value  of  kinematic 
hardening  variable 

Kinematic  hardening  exponent 

Coefficient  of  hardening  recovery 

Hardening  recovery  exponent 

Overstress  parameter 

Strain  rate  sensitivity  parameter 

Saturation  value  of  isotropic 
hardening  variable 

Isotropic  hardening  exponent 


_ Value _ 

24.73  x  103  ksi 
0.336 
30.00  ksi 


350.0 

0.4  x  10-10 


155.0  ksi  sec1'''1 


50.00  ksi 


'  *'4'’ 
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Fig.  13.17  —  Predicted  loading-unloading-reloading  response  at  different  unloading  rates 

by  Bodrter's  theory. 


CONCLUSION 


There  is  no  significant  advantage  to  be  gained  by  abandoning  the  yield  condition-yield  function- 
yield  surface  description  of  non-recoverable  deformation,  whether  it  be  rate  dependent  or  not.  In  a 
continuum  mechanics  description  of  inelastic  behavior  which  views  deformation  on  a  gross 
phenomenological  scale,  the  yield  condition  and  the  yield  function  are  the  means  by  which  this 
manner  of  description  (a)  identifies  that  there  has  been  a  change  in  material  response  due  to  micro- 
structural  changes,  and  (b)  characterizes  that  change. 

The  seeming  simplicity  that  was  to  result  with  not  having  to  consider  a  yield  condition  and  a 
yield  function  and  its  evolution  with  load  history  (hardening),  actually  turns  out  to  provide  no  real 
simplification  at  all.  On  the  contrary,  in  the  absence  of  the  simpler  inherent  theoretical  structure  that 
clearly  links  and  runs  through  the  continually  generalized  line  of  yield  based  viscoplastic  theories  of 
Part  I,  the  proposed  viscoplastic  theories  of  Part  H  resort  to  a  range  of  different,  and  at  times  spuri¬ 
ous,  arguments  that  serve  as  the  foundations  of  each  of  the  different  theories.  These  include:  incor¬ 
porating  microstructural  models  for  polycrystalline  inelastic  deformation  together  with  internal  vari¬ 
ables  into  the  continuum  theory,  where  the  internal  variables  are  loosely  identified  with  the  variables 
of  the  microstructurai  model;  appeal  to  empirical  data  correlation  formulae;  modification  of  descrip¬ 
tive  mechanical  models  for  viscoelastic  behavior;  and  use  of  undefined  ‘hidden'  internal  variables. 

Several  other  viscoplastic  constitutive  theories  have  recently  been  proposed  [67-711  that  were  not 
included  in  this  study.  The  conclusion  drawn  above,  however,  would  not  have  been  altered  had  they 
been  included. 

A  major  difficulty  confronting  viscoplastic  theories  without  the  yield  condition  arises  when 
attempting  to  model  anisotropic  hardening.  As  a  material  hardens,  it  is  actually  continuously  chang¬ 
ing  its  yield  stress  value  as  deformation  proceeds.  The  yield  surface  offers  a  natural  and  simple 
geometrical  representation  of  this  phenomenon,  which  in  turn  provides  the  simplest  means  of  analyti¬ 
cally  describing  or  modeling  the  hardening  process.  This  is  brought  out  by  the  tortuous  and  dubious 
attempt  by  Bodner  and  Stouffer  to  incorporate  so-called  directional  (kinematic)  hardening  into  the 
fabric  of  their  constitutive  theory  This  is  also  shown  explicitly  by  the  difficulty  that  their  theory 
encounters  when  attempting  to  model  uniaxial  lead  cycling.  The  difficulties  experienced  with  uniaxial 
load  cycling  will  probably  become  more  pronounced  for  multiaxial  cyclic  behavior  where  cross¬ 
hardening  effects  enter,  and  for  arbitrary  cyclic  loading. 

The  modified  Chaboche  viscoplastic  theory  appears  to  offer  the  most  promise  for  being  able  to 
model  reasonably  well  the  entire  range  of  inelastic  material  behavior  for  small  deformation. 

Extension  of  vtscoplastic  constitutive  theory  into  the  domain  of  large  deformation  is  possible  in 
principle.  At  the  present  tune,  however,  a  generalization  of  this  kind  that  would  include  anisotropic 
hardening  in  explicit  terms  is  not  as  yet  clearly  understood.  Controversy  (and  thus  confusion)  exists 
over  (a)  how  to  properly  describe  kinematically  large  scale  non-recoverable  deformation,  particularly 
when  there  is  load  cycling,  and  (b)  the  lack  of  frame  indifference  for  the  kinematic  hardening  descrip¬ 
tion  used  for  small  deformation  theory. 

From  the  results  of  the  comparative-qualitative  study  given  in  Part  III,  as  well  as  from  other 
studies,  it  appears  as  if  it  is  now  possible  to  reasonably  model  a  wide  range  of  inelastic  time- 
dependent  material  behavior  characteristics.  Further  development  of  viscoplastic  constitutive  theory 
has  the  potential  for  application  to  the  problem  of  developing  a  general  theory  of  fracture. 
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Appendix 

STRESS,  DEVIATOR  STRESS  AND  STRAIN  MATRICES  FOR  UNIAXIAL  LOADING 


ru  -  a,  T2 2  -  T 33  -  r |2  -  7,3  -  T2 3  -0. 
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a  0  0 
0  0  0 
0  0  0 


(A.  1) 


Correspondingly, 
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from  which 


(A. 2) 
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=  (Tjk  T*,)1/2  =  T2u  +  T2:2  +  T2j3  +  2T2,2  +  2T2,3  +  2T2;3 


(A.  3) 


For  the  elastic  strain:  £«« 


38  2  22  —  pEm  —  Vt' ,  CiV  —  EiV  =  £■»*'  “  0. 


Assuming  inelastic  incompressibility,  TE  =  EM"  +  -22”  ■+"  — 33 ' '  ~  0-  Then  as  -22  + 

E33"  =  -  and  Z12"  =  S33"  (symmetry  of  loading),  -  -22"  =  “33”  = 

—  =  — — e".  Also.  because  of  the  lack  of  shear  stress  components 
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Because  of  the  inelastic  incompressibility 
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